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ABSTRACT 


The high abundances of Complex Organic Molecules (COMs) with respect to 
methanol, the most abundant COM, detected towards low-mass protostars, tend 
to be underpredicted by astrochemical models. This discrepancy might come 
from the large beam of the single-dish telescopes, encompassing several com¬ 
ponents of the studied protostar, commonly used to detect COMs. To address 
this issue, we have carried out multi-line observations of methanol and several 
COMs towards the two low-mass protostars NGC1333-IRAS2A and -IRAS4A 
with the Plateau de Bure interferometer at an angular resolution of 2 arcsec, 
resulting in the hrst multi-line detection of the 0-bearing species glycolaldehyde 
and ethanol and of the N-bearing species ethyl cyanide towards low-mass pro¬ 
tostars other than IRAS 16293. The high number of detected transitions from 
COMs (more than 40 methanol transitions for instance) allowed us to accu¬ 
rately derive the source size of their emission and the COMs colnmn densities. 
The COMs abnndances with respect to methanol derived towards IRAS2A and 
IRAS4A are slightly, but not substantitally, lower than those derived from pre¬ 
vious single-dish observations. The COMs abnndance ratios do not vary signih- 
cantly with the protostellar luminosity, over hve orders of magnitnde, implying 
that low-mass hot corinos are qnite chemically rich as high-mass hot cores. Astro- 
chemical models still underpredict the abundances of key COMs, snch as methyl 
formate or di-methyl ether, snggesting that our understanding of their formation 
remains incomplete. 

Subject headings: astrochemistry - ISM: abundances -ISM: individual objects (NGC 
1333-IRAS2A, NGC 1333-IRAS4A) - ISM: molecules - stars: formation 
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1. Introduction 


The early stages of low-mass star formation are known to be accompanied by the 
increase of the molecular complexity. Most of the lines detected in the sub-millimetric 
spectra of Class 0 protostars are attributed to Complex Or ganic Molecules (COMs, i.e. 


molecules based on carbon chemistry with 6 or more atoms; 


Herbst fc van Dishoeck 


as shown by unbiased spectral surveys of low-mass protostars (see 


Caux et ah 


200911 ■ 


2011 


for the 


spectral survey of the low-mass protostar IRAS 16293-2422 for instance). 


The bright protostars IRAS 16293-2422 and NGC1333-IRAS4A (hereafter IRAS 16293 
and IRAS4A, respectively) have been the hrst two protostars where COMs, such as methyl 


formate, di-methyl ether, fo rmic acid, met 


with single-dish telescopes (jCazaux et ah 


lyl cyanide, or ethyl cyanide , have been detected 


2003 


Bottinelli et ah 


2004aj). The subsequent 


detection of a few COMs towar ds IRAS 16293 with in terfer ometers, providing better 


angular resolutions of rsj 2 ", by 


Bottinelli et ah 


f 2004bll and 


Kuan et ah 


(1200411 conhrmed 


that most of the COMs emission likely comes from the warm inner region of protostellar 
envelopes, called “hot corinos”: the low-mass counterparts of high-mass hot cores. Since 
then, the number of low-mass protostars showing COMs has increa s ed wi th th e detection o 


(2011, 

2014 

)• 

Palau et ah 

(2011 

) and 

Fuente et ah 


Bottinelli et ah 


(I 2 OO 7 I 1 and 


Oberg et ah 


(1201411 also reported the detection of 


several COMs towards four low/intermediate-mass protostars. 


Thanks to a larger number of transitions detected in the broad bands of their receivers, 
single-dish telescopes have been hrst used to derive the column densities of COMs, allowing 
them to constrain their abundances averaged over a relatively large beam of 10-30 ". 

The abundances of COMs, usually compared to that of their probable mother molecules 
(formaldehyde and methanol, see next paragraph), are found to be relatively high: > 10 
% for methyl formate, di-methyl ether, and formic acid, and ~ 1 % for methyl cyanide 
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and ethyl cyanide fsee iBottinelli et al.l 120071: lOberg et al.l 1201 ll. 120141) although the number 
of constrained abundance ratios remains relatively low and the abundances show some 
scattering between the sources and molecules. 

Warm gas phase chemistry triggered by the sublimati on of the main ice co mponents, 


and of methano 


Charnlev et ah 


in p articular, has been hrst invoked by 


Millar et ah 


( 1991 ) and 


(119921) to explain the pre sence of COMs observed towards the high-mass 


hot core Orion KL by 


Blake et al 


(Il987() . However, more recent laboratory experiments 


and theoretical calculations have contradicted several key assumptions made in the gas 
phase models: dissociative recombination of large ions do not lead pre dominantly to the 


formation of COMs b ut rather to their fragmentation into small pieces (iGeppert et ah 


Hamberg et ah 


2006 


20101 ) whilst ion-molecule reactions have bee n found to be no t sufficiently 


efficient to produce the observed amount of methyl formate (iHorn et ah 


2004 ). The current 


scenario of COMs formation is now based on the recombination of radicals at the surface 


of interstellar grains during the warm-up phase (30 K < T < 10 0 


envelopes surrounding Class 0 protostars (jCarrod fc Herbst 


2006 


K) occurring in t 


Carrod et ah 


20081). 


le 


In these models, the radicals are generated by the UV photodissocation of the main ice 
components or they have survived to the incomplete hydrogenation process of C O leading 


to CH 3 OH during the ice formation (jCarrod fc Herbst 


2006 : 


Taqnet et al. 


2012 ). Since 


the 90s, laboratory experiments have shown that UV irradiation of interstellar ice analogs 
containing methanol, formaldehyde, and am monia can lead to the formation of a pletho ra 


of complex molecules and even amino-acids ( 

All 

amanc 

ola et al. 1 

988; 

Cerakines et al. 

1996 ; 

Hudson & Moore 

2000 ; 

Munoz Caro et al. 

2002 


Oberg et al 

200 < 

3). However, the 


quantitative efficiency of the COMs formation in ices and their actual chemical pathways 
are still highly uncertain. 


Although the current models produce a large set of COMs in significant quantities, they 
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are not able to explain the very high abnndance ratio (> 10 %) with respect to methano l 


seen for a few COMs, snch as methyl formate or di-methyl ether (see 


Tag net et af 


2012 


for a discnssion of this problem). The discrepancy between observations and models conld 
be dne to the large beams of single-dish telescopes nsed to derive the abnndance ratios of 
COM s. Typical single-dish beams of 10 " are mnch larger than the size of hot corinos 


0.5 


Maret et ah 


2005 


Manrv et ah 


20141) and encompass the cold envelopes a nd possible 


ontflows driven by the c e ntral protostars, w here COMs have also been detected flArce et ah 


2008 

Obere: et al. 

201ol 

Jaber et al. 

2014) 


in cold and qniescent cores by 


20141). In addit ion, COMs have also 


oeen found 


Racmann et al. 


1120121) and 


Cernicharo et al 


(2012) who 


claimed that quiescent cold ga s phase chem i stry ca n produce COMs but in lower quantities. 


However, new observations by 


Vastel et ahl (120141) rather suggest that the emission from 


COMs observed in another pre-stellar core originates in an outer ring, so that the previous 
conclusions may need some cautions. 

To better constrain the abundances of COMs originating from the hot corinos 
surrounding the low-mass protostars, a large number of transitions of methanol and COMs 
need to be observed with interferometers, providing angular resolutions of ~ 1-2 ". The 
emission originating from hot corinos can be distinguished from other components of 
the envelope, allowing us to directly derive the abundance of COMs in the hot corinos. 
Moreover, CH.^OH emission is likely optically thick towards the continuum peak of 


protostars fIZapata et al. 


20131) . observations of its optically thin isotopologue ^^CHsOH 


are therefore required to derive an accurate estimate of the methanol column density. 
Alt hong 
corinos 


1 several publications 

lave reported the interferometric d( 

Tection of COMs in 1 

Bottinelli et al. 

2004b: 

Kuan et al. 

2004 

Joreensen et al. 

2005: 

Bisschop et al. 


J 


2008 


Jorgensen et al 


2011 


Persson et al 


201 


Jorgensen et al. 


2012 


Manrv et al 


20141 ). 


to our knowledge, none of them led to an accurate and simultaneous estimation of the 
column densities of methanol and COMs. 
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In this work, we present multi-line observations of methanol (^^CHsOH, and ^^CHsOH) 
as well as several COMs (methyl formate, di-methyl ether, ethanol, glycolaldehyde methyl 
cyanide, and ethyl cyanide) performed with the Plateau de Bure interferometer towards 
the two low-mass protostars IRAS2A and IRAS4A located in the NGC1333 star-forming 
region. Although the angular resolution does not allow us to spatially resolve the emission 
of COMs, it is sufficiently high to distinguish the emission from the hot corinos from other 
components of protostellar envelopes. The paper is structured as follows: section 2 describes 
the observational strategy, section 3 presents the continuum maps as well as the spectra and 
the maps of molecular transitions, section 4 explains the adopted methodology to derive the 
abundances of COMs, section 5 discusses the results; and section 6 summarizes this work 
with the conclusions. 


2. Observations 

The two low-mass Class 0 protostars IRAS2A and IRAS4A were observed with the 
IRAM Plateau de Bure Interferometer (PdBi) on 2010 July 20, July 21, August 1, August 
3, November 24, and 2011 March 10 in the C and D conhgurations of the array. Due to 
the proximity to each other, the two sources were observed in the same track. Phase and 
amplitude were calibrated by performing regular observations of the nearby point sources 
3C454.3, 3C84, and 0333-1-321. The amplitude calibration uncertainty is estimated to 
be ~20%. The WIDEN backends have been used at 143.4 and 165.2 CHz, providing a 
bandwidth of 3.6 CHz each with a spectral resolution of 1.95 MHz (~ 3.5 — 4 km s“^). 
High-resolution narrow band backends focused on two CH 3 OH lines and 12 HCOOCH 3 lines 
have also been used. They provide a bandwidth of 80 MHz with a spectral resolution of 
0.04 MHz (0.08 km s“^). Due to the low signal-to-noise ratio (S/N) obtained for the methyl 
formate lines at high spectral resolution, we decreased the spectral resolution to 0.4 MHz 
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(0.8 km s“^) to obtain a S/N ratio higher than 3. The data calibration and imaging were 
performed nsing the CLIC and MAPPING packages of the GILDAS software 0. Gontinnnm 
images were prodnced by averaging line-free channels in the WIDEN correlator before 
the Fonrier transformation of the data. The coordinates of the sonrce and the size of the 
synthesized beams are reported in Table 1. 


3. Results 


3.1. Continuum Maps 


Figured] shows the maps of the continuum emission of IRAS2A and IRAS4A at 143 and 
165 GHz obtained after natural weighted cleaning. Parameters of the continuum emission 
(integrated flux and deconvolved FWHM size), obtained from elliptical Gaussian hts in the 
{u, v) plane, are given in Table dl For the two settings, the FWHM size of the continuum 
emission is slightly smaller than the size of the synthesized beam, the continuum emission is 


consequently not resolved. 
" separation (iLooney et ah 


n pa rticular. IRAS4A is known to be a binary system with a 1.8 


20001), as depicted by the two red crosses in Fig. dlthat indicate 


the positions of IRAS4A-SE and -NW. Although the continuum emission of IRAS4A is 
peaked at the southeast (SE) position rather than at the northwest (NW) position for the 
two settings, we cannot resolve the two sources. 

We estimated the mass of the envelope and the H 2 column density from the continuum 
fluxes and the sizes derived from elliptical Gaussian hts listed in Table dl Assuming an 
optically thin dust emission, the mass M of the envelope is given by 


M = 




K^By{Td)Rd 


( 1 ) 


^http://www.iram.fr/IRAMFR/GILDAS 










Table 1: Properties of NGC 1333 IRAS2A and IRAS4A. 


Source 

1RAS2A 

1RAS4A 

R. A. (J2000) 

03:28:55:57 

SE: 03:29:10.52 



NW: 03:29:10.42 

Bed. (J2000) 

31:14:37:22 

SE: 31:13:31.06 



NW: 31:13:32.04 

d (pc) “ 

235 

235 

Ulsr (km/s) 

+7.7 

+7.2 

Lboi (Lq) ^ 

36 

9.1 

Menv (Mq) “ 

5.1 

5.6 

Frequency = 143 GHz 

Beam size (") 

2.2x1.7 

2.1x1.7 

Beam PA (°) 

25 

25 

rms(WideX) ^ 

2.57 

3.34 

rms(Cont.) ® 

1.56 

10.8 

Flux (Jy) 

0.13 

1.1 

Size (") ^ 

1.7x1.7 

2.1x1.7 

PA (°) t 

+51 

+25 

M (Mq) 9 

0.4 

3.8 

N{R2) (cm-2) a 

5.0 X 10^4 

3.7 X 10^5 

Frequency = 165 GHz 

Beam size (") 

2.3x1.7 

2.4X1.8 

Beam PA (°) 

no 

115 

rms(WideX) ^ 

3.50 

4.02 

rms(Gont.) ® 

1.84 

10.8 

Flux (Jy) t 

0.19 

1.6 

Size (") f 

1.9x1.7 

2.5X1.6 

PA (°) t 

-67 

-44 

M (Mq) a 

0.4 

3.5 

N{ll2) (cm-2) 3 

4.5 X 10^4 

3.1 X 10^5 


Note. — a: 


Hirota et al 


iaOOSl); b: 


Karska et al 


(l2013h : 


Kristensen et al. 


( 2012 1: d: Units of 


mjy/beam/channel for a channel width of 1.95 MHz. e: Units of mJy/beam. /: Continuum integrated 


fluxes and sizes were obtained from elliptical Gaussian fits in the (u, v) plane (i.e., deconvolved FWHM size); 
o: The envelooe mass and avera.ffed column densitv were derived from the continnnm fluxes obta.ined within 
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Fig. 1.— Continuum maps at 145 and 165 GHz of IRAS2A (rms of 1.56 and 1.84 mJy/beam, 
respectively) and IRAS4A (rms of 10.8 and 10.8 mJy/beam, respectively) . The positions of 
the sources IRAS2A, IRAS4A-NW, and IRAS4A-SE are marked by a red plus sign. Contour 
levels are in steps of 3cr. The white ellipse represents the half-power beamwidth of the 
synthesized beam. 


where S,^ is the continuum flux integrated over the gaus sian ellipse and list ed in Table [H 


d is the distance to the two low-mass protostars (235pc; 


Hirota et ah 


2008), By{Td) is the 


Planck black body function for a temperature assumed to be 30 K, Rd is the dust-to-gas 
m ass ratio equal to 0.01. K„ is the opacity per dust mass taken from column 5 in Table 5 


of lOssenkopf fc Hennind (jl994f ) (corresponding to grains showing a MRN size distribution 
covered by thin ice mantles at uh = 10® cm“^) and extrapolated to 145 GHz = 0.38 
g/cm^) and 165 GHz = 0.48 g/cm^). The column density of H 2 averaged over the 
gaussian ellipse can be deduced from M following this formula 


iV(H2) 


M 

lirriYiVtd? 


( 2 ) 


where M is the envelope mass, /r = 2.38 is the mean molecular mass in units of hydrogen 
atom masses, mn is the hydrogen atom mass, and hi the solid angle subtended by the 
gaussian ellipse. The envelope mass M and H 2 column density A^(H 2 ) derived for the two 
frequencies are listed in Table [H 
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3.2. Spectra 


For the two sources and for all the settings, we obtained the spectral cubes by 
subtracting the continuum visibilities to the whole (line+continuum) datacube. For 
IRAS4A, the baseline has been flattened by importing the data cubes into CLASS and 
subtracting a polynomial function to each individual spectrum. The 3.6 GHz wide spectra 
obtained at the coordinates of IRAS2A and IRAS4A-NW with the two WideX backends are 
presented in Figure [2J We also present the narrow band spectra of CH 3 OH and HCOOCH 3 
in Figures |3] and m The la rms noise in the line-free channels of all spectra are given in 
the caption of Fig. [2] and [3] and HI The two sources are chemically rich since the WideX 
spectra towards IRAS2A and IRAS4A-NW display ~ 200 and ~ 170 lines detected with a 
signal-to-noise ratio higher than 3, resulting in a l ine de nsity of 28 and 23 detected lines per 


GHz, respectively. For comparison. 


Maury et al.i (120141) detected 86 lines above the 3-sigma 


level between 216.9 and 220.5 GHz towards IRAS2A with the PdBi, resulting in a similar 
line density of 24 detected lines per GHz. 


Th e line identihcation was carried out using the JPL fjPickett et al 


19981) and the 


GDMS fiMiiller et ahl 120051) spectroscopic catalogues. Line identideations were performed 


by eye, by taking into account the upper energy level, the line strength, and the velocity of 
each transition. The detected molecules are listed in Table 2 with the number of detected 
transitions, the energy range of their upper energy levels, and the spectroscopic reference. 
We detected about 35 transitions for the main isotopologue of methanol with upper energy 
levels up to ~ 1020 K and about 13 lines for its isotopologue ^^GH 30 H. In addition to 
methanol, we report the hrst multi-line detection of glycol aldehyde HGOGH 2 OH, ethyl 


cyamd 

16293 


e G 2 H 5 GN, and ethanol G 2 H 5 OH towards low-mass protostars other than IRAS 




We also detected several transitions originating from methyl formate HGOOGH 3 , 


^ The detection of glycolaldehyde was reported in IRAS2A by Goutens et al. (2015) while 













Flux [Jy/beam] Flux [Jy/beam] 
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163.5 164 164.5 165 165.5 166 166.5 167 

Rest Frequency [GHz] 

Fig. 2.— PdBi continuum-subtracted spectra of the WideX backends around 143 and 165 
GHz towards the peak position of IRAS2A and IRAS4A-NW. The rms noise levels are 2.57 
and 3.50 mJy/beam/channel at 143 GHz and 165 GHz towards IRAS2A and 3.34 and 4.02 
mJy/beam/channel at 143.5 GHz and 165 GHz towards IRAS4A, respectively. 
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Fig. 3.— PdBi spectra of methanol obtained with the narrow band correlators towards 
the peak position of IRAS2A and IRAS4A-NW The rms noise levels are 7.80 and 7.86 
mJy/beam/channel towards IRAS2A and IRAS4A, respectively. 



Fig. 4.— PdBi spectra of methyl formate obtained with the narrow band correlators to¬ 
wards the peak position of IRAS2A and IRAS4A-NW. The rms noises are about 7-10 
mJy/beam/channel towards IRAS2A and IRAS4A. 
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di-methyl ether CH 3 OCH 3 , and methyl cyanide CH 3 CN towards the two sonrces. Dne to 
the low spectral resolntion, we carefnlly checked that the detected lines do not snffer from 
any blending with other transitions from similar or other molecnles by nsing the JPL and 
CDMS databases bnt also with the Splatalogne database 1^. In total, abont 70 % of the 
lines detected at a 3 sigma level in the two 3.6 GHz WideX correlators are attribnted to the 
complex organics listed in Table 2. Other lines are attribnted to the denterated methanol 
isotopolognes CH 2 DOH, CH 3 OD. CHD^OH, ot her denterated molecnles snch as HDO 


(stndied in a previons work; 


Tag net et ah 


2013 1. DCN, NH 2 D, and DCO"^ (in absorption) 


and the snlphnr-bearing species SO 2 , and C^^S. The analysis of the denterated methanol 
transitions will be pnblished in a separate article. ~20 lines are nnidentified. 

We estimated the Fnll-Width-at-Half-Maximnm (FWHM) of the lines detected with 
the narrow band correlators giving a spectral resolntion of 0.4 MHz, and with the WideX 
correlator providing a spectral resolntion of 1.95 MHz, throngh a ganssian £t of the spectra 
obtained at the coordinates of IRAS2A and IRAS4A-NW. Table [3] lists the linewidths of 
the transitions detected with the two correlators. The nncertainties in the linewidths are 
dne to the statistical errors from the ganssian fit and to the nncertainty from the low 
spectral resolntion. In this work, the low spectral resolntion dominates the nncertainty 
of the linewidths. Tables [S] to [H] of the Appendix list the properties of all detected 
transitions along with the FWHM linewidths. For the two sonrces, the FWHM of the 
CH 3 OH transitions detected at high spectral resolntion are abont 3 km/s whilst the widths 
of the HCOOCH 3 lines vary between 0.9 and 2.8 km/s when the lines of the -E and 


the -A states do not overlap. The linewidths derived with the Pc 
linewidths of other CH 3 OH and HCOOCH 3 transitions derived by 


Hi are similar to the 


Maret et ah 


(1200511 and 


this article was in the review process, 
^http: / / WWW. cv.nrao.edn / php / splat / 
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Table 2: List of molecules detected towards IRAS2A and IRAS4A. 



1RAS2A 


IRAS4A 

Ref. 

Molecule 

-^lilies 

Eup 

(K) 

-^lines 

(K) 

CH 3 OH 

34 

14 - 1022 

35 

14 - 1022 

1 

HCOOCH 3 

20 

43 - 248 

20 

43 - 237 

2 

CH 2 DOH 

13 

33 - 230 

13 

33 - 230 

3 

13 CH 3 OH 

13 

14 - 222 

12 

14 - 222 

4 

CH 3 OCH 3 

8 

11 - 314 

7 

11 - 314 

5 

C 2 H 5 OH 

8 

37 - 216 

7 

37 - 216 

6 

CH 3 CN 

7 

40 - 390 

7 

40 - 390 

7 

CHD 2 OH 

6 

20 - 67 

6 

20 - 67 

8 

HCOCH 2 OH 

4 

53 - 68 

7 

53 - 177 

9 

SO 2 

4 

24 - 102 

4 

24 - 102 

10 

C 2 H 5 CN 

4 

63 - 130 

3 

63 - 130 

11 

CH 3 OD 

1 

40 

1 

40 

12 

H 2 I 3 CO 

1 

10 

1 

10 

13 

H 2 C 1®0 

1 

22 

1 

22 

14 

HC 3 N 

1 

75 

1 

75 

15 

NH 2 CHO 

1 

30 

1 

30 

16 

CH 2 CO 

1 

41 

1 

41 

17 

C34S 

1 

14 

1 

14 

18 

DC 3 N 

1 

62 

- 

- 

15 

HDO 

1 

319 

1 

319 

19 

DCN 

1 

10 

1 

10 

20 

NH 2 D 

1 

183 

1 

183 

21 

D 2 CO 

1 

21 

1 

21 

14 


Note. — 

1: 

Xu et al. 

5: 

Lovas et al. 

(19791 

; 6 

Carrol et al. 

(2 

OIC 

); 10: 

Muller et al. 

(2 

OOf 

); 14: 

Fabricant et al 

(1977 

): 

Cohen & Picke 

tt ( 

198 

2); 


(I2nn8h: 2: 


Alekseev et al 


Ilvus 


Endres et al. 


linetal 


(2009); 7: 


Dangoisse et a' 


Bogey et al 


(2009): 3: 


Pearson et al. 


( 201211 : 4: 


Cazzoli and Puzzarini (1200 


1996): 11: 


(ll97Sh : 15: 


Pular2ania^t^ 


(II 982 I 1 : 19: 


La^ert^&Lq^ 


Messer et al 




(199611: 12: 


Xu fc Lovas ( 19971 1: 


Parise et al 


Anderson et al 


197811 : 16 : 


(1198411 : 20: 


Johnson et al 


De Lucia fc Gordv 


( 1988 1: 13 


( 2 OO 2 II : 9 


( 1972 1: 17 


1969 1: 21 
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Bottinelli et al. 


fj2004al) . Due to the low spectral resolution, the linewidths derived from the 


WideX correlators result from the convolution of the intrinsic linewidths with the spectral 
resolution of 1.95 MHz. The top panel of Figure [5] shows the deconvolved FWHM of 
CH 3 OH, ^^CHsOH, HCOOCH 3 , and CH 3 CN, showing a high number of transitions detected 
at a high signal-to-noise ratio, as a function of the energy of the upper level towards the 
two sources. For this purpose, we excluded several CH 3 OH and HCOOCH 3 transitions 
that are blended with each other. No clear trend can be deduced for the two sources. The 
fluctuation of the deconvolved FWHM linewidths between 2 and 8 km/s seems to be due 
to their high uncertainties. The bottom panel of Fig. |5] compares the linewidths deduced 
from the narrow band correlators giving a high spectral resolution of 0.4 MHz with the 
deconvolved linewidths from the WideX spectra for the lines detected with the narrow band 
correlators. For all the lines, the FWHM linewidth deconvolved from the WideX spectra 
is higher than the FWHM linewidth deduced from the narrow band correlator but the 
differences remain within the uncertainties. 


3.3. Line Maps 

For all the transitions, the interferometric maps of the IRAS2A and IRAS4A protostars 
have been obtained by integrating the flux over 14 sR=t^^ where Wsr is the system velocity 
of the source and AV = 3 km/s following the FWHM linewidths of the CH 3 OH transitions 
listed in Table [3l In practice, dne to the low resolntion of the WideX backends, the line 
emission is integrated over 3 channels. Figures |6] and [7] show a compilation of the integrated 
line maps towards IRAS2A and IRAS4A obtained after natnral weighted cleaning. For 
species where several transitions were detected, two maps showing a low-energy and a 
high-energy transition are presented. Tables |8] to [16] of the Appendix list the properties 
of all detected transitions along with their FWHM (fnll-width-at-half-maximnm) sizes 
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Table 3: Linewidths of the transitions detected in the narrow band correlators. 


Molecule 

Frequency 

(GHz) 

E 

■^up 

(K) 

IRAS2A 


IRAS4A 

dVfi 

(km/s) 

dVw 

(km/s) 

dVn 

(km/s) 

dVw 

(km/s) 

CH 3 OH 

143.86580 

28.3 

2.9 ± 0.8 

6.7 ± 4.1 

3.2 ± 

0.8 

7.5 ± 4.1 

CH 3 OH 

143.16952 

113 

3.2 ± 0.8 

6.1 ± 4.1 

2.8 ± 

0.8 

5.8 ± 4.1 

HCOOCH 3 

143.23420 

47.3 

2.3 ± 0.8 

5.8 ± 4.1 

1.6 ± 

0.8 

4.1 ± 4.1 

HCOOCH 3 

143.24051 

47.3 

2.8 ± 0.8 

5.8 ± 4.1 

1.9 ± 

0.8 

4.1 ± 4.1 

HCOOCH 3 

163.82968 

62.5 

2.4 ± 0.7 

6.7 ± 3.6 

2.0 ± 

0.7 

4.7 ± 3.6 

HCOOCH 3 

163.83553 

62.5 

1.9 ± 0.7 

3.6 ± 3.6 

2.4 ± 

0.7 

4.4 ± 3.6 

HCOOCH 3 

163.96039 

64.5 

5.3 ± 0.7 

6.9 ± 3.6 

4.1 ± 

0.7 

5.9 ± 3.6 

HCOOCH 3 

163.96188 

64.5 






HCOOCH 3 

163.98746 

64.5 

5.2 ± 0.7 

6.9 ± 3.6 

0.9 ± 

0.7 

6.0 ± 3.6 

HCOOCH 3 

163.98891 

64.5 



2.2 ± 

0.7 


HCOOCH 3 

164.20598 

64.9 


5.9 ± 3.6 

0 

bo 

0.7 

6.6 ± 3.6 

HCOOCH 3 

164.22382 

64.9 

2.1 ± 0.7 

6.5 ± 3.6 

2.3 ± 

0.7 

5.2 ± 3.6 


Note. — dVn and dVw are the non-deconvolved linewidths deduced from the spectra obtained with the 


narrow band and WideX correlators, respectively. 
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of emission and their Position Angle (PA) derived from the modelling of the visibilities 
assuming elliptical gaussians, circular gaussians or point sources if gaussian hts were not 
possible. 


For the two sources, the emission of most molecules is limited to the inner regions 
near the protostars. In IRAS4A, the compact emission of all COMs transitions originates 
from the -NW source although the SE protostar is brighter in the continuum. In contrast, 
no molecular lines seem to originate from IRAS4A-SE, as also obs erved in previ o us 


interferometric observations of II 2 O and other complex organics by 


Persson et ah 


(I2OI2I). 


The low-energy transitions of CII 3 OII (the transition at ~ 145.094 GHz including several 
non-resolved transitions towards IRAS2A not shown in this work and all the transitions 
with Eup < 120 K towards IRAS4A) show an extended emission consistent with the position 
of molecular outflows. For IRAS2A, the interferometric map of the CH 3 OH transitions 
at ~ 145.094 GHz displays a slightly red-shifted emission located 7 " to the north from 
IRAS2A, co nsistent with the direc tion of an outflow previously detected at small scales 
with GO by 


Jorgensen et ah 


(120071) . For IRAS4A, the emission of the low-energy GH 3 OH 


lines extends in the bipolar outflow along a N-S direction and seems to peak at ~ 10 


south from the protostars, consistent with the south lobe detected in SiO by 


Ghoi et ah 


(120051 ). From Fig. [TJ it is clear that the inner hot corino does not dominate the 


lux of 


the weakly excited transitions of methanol in IRAS4A (see also 


Maret et af 


20051) . The 


emission of the HG 3 N transition is also spatially resolved and shows an elongation of its 
emission towards the SW direction for IRAS2A and towards the N direction for IRAS4A. 
The spatial distribution and the kinematics of the outflow driven by IRAS2A and IRAS4A 
will be analysed in a future work. 


The emission originating from all GOMs transitions but some GH 3 OH and HG 3 N lines 
is not spatially resolved by the array since their FWHM emission size deduced from the £t 
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to the visibilities is lower than the synthesized beam of the interferometer. Conseqnently, 
the sonrce sizes presented in this work can only be nsed as npper limits. These observations 
are qnalitatively consistent with previons models and observations snggesting that methanol 
and COMs mostly come from the inner hot corino. In this region, these molecnles show 
a jnmp of their abnndan ce when the t empe ratnre is higher than the temperatnre of ice 


snblimation (T ~ 100 K). 


Maret et ah 


(I 2 OO 4 J) estimated a size for the hot corinos of IRAS2A 


and IRAS4A of 0.45 " by reprodncing the formaldehyde emission observed with single¬ 
dish telescopes with an abnndance jnmp of two orders of magnitnde at r ~ 5 0 AU from 
the ce ntral protostar. The Inminosities assnmed for IRAS2A and IRAS4A by 


Maret et al 


(j2004j) are lower by a factor of 2.25 and 1.5, respectively, than the Inminosities assnmed in 


this work. Assnming that the temperatnre profile is governed by the Stefan-Boltzmanns 
law implies a difference in the temperatnre of a factor of 1.2 and 1.1, respectively. The 
assnmption of higher In minosities for the tw o sonrces wonld therefore increase the size of 


the corino of a few AU. 


Maury et al. 


( 20141) estimated a FWHM size of 0.4-0.9" for the hot 


corino of IRAS2A through the use of the more extended A configuration of the PdBi array. 


4. Chemical Abundances in the Hot Corinos 
4.1. Rotational Diagrams 

First estimates of the excitation temperatures and the column densities of observed 
molecules in the hot corinos of IRAS2A and IRAS4A have been obtained from the rotational 
diagram (RD) analysis by assuming optically thin emission and a LTE population of the 
levels. Since the emission of most transitions is not spatially resolved with the PdBi, 
we measured the flux of all transitions originating from a circular mask with a diameter 
equal to the major axis of the synthesized beam size of the telescope (~2.1-2.6"). Tables 
[8] to [16] of the Appendix list the measured flux for all transitions and in the two sources. 
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Astrochemical/dynamical models predict low abundances of COMs in the the dense regions 
of protostellar envelopes where the dust temperature is lower than the temperature of ice 
sublimation Tgv 100 K) due to the efficient depletion, while the abundan ce profiles 


show a strong jump once T = Tev (see 


Aikawa et ah 


2008 


Tag net et ah 


20i 4 induced by 


thermal evaporation in the hot corino. We therefore assumed that all the flux measured 
in the ~ 2 " mask comes from the hot corino region. Although the emission sizes of the 


COMs transitions are not necessarily similar (in IRAS 16293 for exa mple, methy 


has a source size of 4 " whereas formamide has a source size of 2 " 


Jaber et al. 


formate 


201411 


all the COMs observed towards the two sources (iMaret et ah 


assumed a ho 

; corino size 6^ of 0.5 " for 

(Maret et al. 

2004 

). The linewidth at 


FWHM AV is f ixed to 3 krn / s, wh ich r epresents an average va lue of Table [3] and previous 


observations by 


Maret et ah 


fcoOhl l and 


Bottinelli et ah 


(120071) . For molecules with only 


a few lines detected within a narrow range of upper energy levels (glycolaldehyde, or 
ethyl cyanide) or for molecules showing only one detection, we assumed two values for the 
rotational temperature Trot: Trot = Trot(CH 30 H) and Trot = Tevap(ice) = 100 K. As already 


shown in several published observat ional works 


study i ng the emission of methanol towarc 


2013 


Zapata et al. 


(Parise et al. 

2006; 

Bisschop et al. 

2007 


Isokoski et al 


2013fl . it is likely that low-energy transitions of CH 3 OH are optically 


thick towards the center of protostellar envelopes, giving rise to an underestimation of their 
population in the rotational diagrams. Consequently, we excluded the CH 3 OH transitions 
with Tup < 200 K from our RD analysis. The rotational temperatures and total column 
densities of all species derived towards the two sources are summarized in Table lH 


For most molecules, observational data can be reasonably well fitted by a straight line 
with some scattering, likely due to opacity or non-LTE effects. For most of species, the 
column densities derived towards IRAS2A and IRAS4A are very similar. Column densities 
of CH 3 OH are about 6 — 12 x 10^^ cm“^ whilst the column densities of complex organics 
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range between 6 x 10^^ and 6 x 10^® cm”^. The rotational temperatnres derived for most 
COMs are generally higher than 100 K in IRAS2A and in IRAS4A. Althongh they do not 
not necessarily reflect the kinetic temperatnres, the high rotational temperatnres fonnd 
in the inner protostellar envelope are in good agreement with the kinetic temperatnres 
expected in hot corinos (T > 100 K). The CH 3 CN popnlation distribntion can be nsed to 
estimate the kinetic temperatnre of the warm inner envelope becanse the CH 3 CN popnlation 
distribntion over the different K ladders observed f or J = 9 can only be modihed by 


collisions as radiative transitions are prohibited (see 


Wang et ah 


2010 


, for a more detailed 


discnssion of the CH 3 CN popnlation in hot cores). If the energy states are thermalized 
and the transitions are optically thin, the kinetic temperatnres within the hot corinos of 
IRAS2A and IRAS4A wonld be close to 290 and 360 K, respectively. As we will see in the 
next section, the temperatnres are probably overestimated since the CH 3 CN transitions are 
likely optically thick. 

We derived a abnndance ratio of CH 3 OH of 26 and 12 in IRAS2A and 

IRAS4A, respectively. The two yalnes are lower than the abnndance ratio of 


70 expected in the local ISM (IBoogert et ah 


2002 


Milam et ah 


2005 : 


Wirstrom et ah 


20 111 ) by a factor 2.7 and 5.8, respectively. We verihed that nsing the same rotational 
temperatnre for CH 3 OH and ^^CIl 30 II only modihes the abnndance ratios by a few 

percents at most. The low abnndance ratio might be dne to the different ranges 

of excitation of the observed transitions nsed to derive the colnmn densities. We nsed 
excited transitions (with E^p > 200 K) to derive W(^^CIl 30 H) whereas only weakly excited 
^^CHsOII transitions (with E^p < 225 K) have been detected. The low ratio measnred 
in the two protostars snggests that excited transitions of CII 3 OH {E-^p > 200 K) are also 
optically thick. An overestimation of the ^^CIl 30 II colnmn density by the RD best £t, dne 
to the large nncertainty on the flnxes, is also possible. An analysis taking the opacities into 
acconnt is therefore reqnired to clarify this issne. 
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Table 4: Results from the rotational diagram analysis. 


Molecule 

Nhc “ 

'Trot 

X " 

-^meth 

NhASD) 

V„eth(SD) 


(cm“^) 

(K) 


(%) 

(cm“®) 

(%) 

IRAS2A 

CH3OH ^ 

(1.2 ± 0.4) X 10^® 

179 ± 62 

(2.5 ± 0.9) X 10“’’ 


1.4 X lO^’’ 


^^CHaOH 

(4.8 ± 1.3) X 10^® 

164 ± 43 

(9.6 ± 2.5) X 10“'* 




HCOOCH3 

(6.4 ± 1.9) X 10^® 

200 ± 61 

(1.3 ± 0.4) X 10“® 

1.9 ± 0.8 

< 1.2 X 10^'^ 

< 85 

CH3CN 

(1.0 ± 0.2) X 10^® 

289 ± 63 

(2.0 ± 0.4) X 10“® 

0.30 ± 0.10 

1.5 X 10’-'^ 

1 

CH3OCH3 

(4.1 ± 1.6) X lO^*^ 

154 ± 62 

(8.2 ± 3.3) X lO”'* 

1.2 ± 0.6 

< 7.2 X 10^® 

< 53 

C2H5OH 

(5.1 ± 2.2) X 10^*^ 

325 ± 140 

(1.0 ± 0.4) X 10“* 

1.5 ± 0.8 

/ 

/ 

HCOCH2OH 

7.8 X 10^® 

179 f 

1.6 X 10“® 

0.23 ± 0.06 

/ 

/ 


2.5 X 10^® 

100 f 

5.0 X 10“^® 

0.074 ± 0.019 



C2H5CN 

1.2 X 10^® 

179 f 

2.4 X 10“^® 

0.036 ± 0.010 

< 1.7 X 10^® 

< 13 


6.9 X 10^^ 

100 f 

1.4 X 10“^® 

0.021 ± 0.005 



HC3N 

7.0 X 

179 f 

1.4 X 10“1® 

0.021 ± 0.005 

/ 

/ 


7.1 X 10^“^ 

100 f 

1.4 X 10“^® 

0.021 ± 0.006 



H2^^C0 

6.6 X 10^® 

179 f 

1.3 X 10“® 

0.20 ± 0.05 

/ 

/ 


2.1 X 10^® 

100 f 

4.3 X 10“^® 

0.063 ± 0.017 



NH2CHO 

1.2 X 10^® 

179 f 

2.3 X 10“® 

0.35 ± 0.09 

/ 

/ 


4.3 X 10^^ 

100 f 

8.7 X 10“^® 

0.13 ± 0.03 



CH2CO 

7.0 X 10^® 

179 f 

1.4 X 10“® 

0.21 ± 0.05 

/ 

/ 


2.6 X 10^® 

100 f 

5.2 X 10“^® 

0.077 ± 0.020 



IRAS4A 

CH3OH ® 

(6.3 ± 3.1) X 10^’’ 

300 ± 151 

(1.7 ± 0.9) X 10“* 


2.0 X 10^’’ 


^^CHaOH 

(5.1 ± 1.5) X 10^® 

197 ± 56 

(1.4 ± 0.4) X 10“® 




HCOOCH3 

(5.2 ± 3.3) X 10^® 

141 ± 90 

(1.4 ± 0.9) X 10“® 

1.5 ± 1.0 

1.0 X lO'^’’ 

52 

CH3CN 

(6.5 ± 2.9) X 10^® 

360 ± 162 

(1.8 ± 0.8) X 10“^® 

0.18 ± 0.10 

2.6 X lO'^'' 

1 

CH3OCH3 

(3.1 ± 1.0) X 10^® 

86 ± 27 

(8.5 ± 2.6) X 10“1® 

0.87 ± 0.37 

< 4.5 X 10^® 

< 22 

C2H5OH 

(4.4 ± 1.4) X 10^® 

221 ± 69 

(1.2 ± 3.7) X 10“® 

1.2 ± 0.5 

/ 

/ 

HCOCH2OH 

(8.9 ± 3.4) X 10^^ 

124 ± 48 

(2.4 ± 0.9) X 10“1® 

0.25 ± 0.12 

/ 

/ 

C2H5CN 

2.3 X 10^® 

300 f 

6.2 X 10“11 

0.064 ± 0.018 

1.9 X 10^®^ 

< 0.92 


8.2 X 10^^ 

100 f 

2.2 X 10“^^ 

0.023 ± 0.007 



HC3N 

6.7 X 101“^ 

300 f 

1.8 X 10“11 

0.019 ± 0.005 

/ 

/ 


6.8 X 10^“^ 

100 f 

1.8 X 10“^^ 

0.019 ± 0.005 



Ha^^CO 

5.3 X 10^® 

300 f 

1.4 X 10”^® 

0.15 ± 0.04 

/ 

/ 


1.1 X 10^® 

100 f 

3.1 X 10“^^ 

0.032 ± 0.009 



NH2CHO 

8.5 X 10^® 

300 f 

2.3 X 10“1® 

0.24 ± 0.07 

/ 

/ 


2.1 X 10^® 

100 f 

5.7 X 10“^^ 

0.059 ± 0.017 



CH2CO 

1.6 X 10^® 

300 f 

4.2 X 10“^® 

0.43 ± 0.12 

/ 

/ 


4.0 X 10^® 

100 f 

1.1 X 10“^® 

0.11 ± 0.03 




Note. — 

a: Column densities averaged over a source size of 0.5 "(see text), b: The abundances relative to H 2 
are obtained from A^(H 2 ) derived at 145 GHz in Table [1] assuming an homogeneous H 2 column density 
within the beam, c: The abundance ratios are relative to the ^^CHaOH column density multiplied by 70 
(see text), d: Colu mn densities derive d from previous single-dish observations an d scale d to a source size 


of 0.5 ". CHaOH: 


IRAS4A: 


Maret et al 


Bottinelli et al 


(l2005[l . other COMs in IRAS2A: 


Bottinelli et al 


( 2007 ). other COMs in 


(|2004al) . e: N^c and T^ot have been derived from the rotation diagram neglecting 


the transitions with E^p < 105 K (see text), f: When they could not be derived from the rotation diagram, 
the rotational temperatures have been assumed to be equal to the temperature of evaporation of water ice 
and the rotational temperature of CH 3 OH. 
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4.2. Population Diagrams 


We used the so-called Population D iagram (PD) analysis following the method 


described by 


Goldsmith fc Langerl fll999l) to investigate the effect of optical depth on the 


column densities of each level. We applied the PD analysis where four or more transitions 
were detected for each species. Briefly, the PD analysis includes the influence of optical 
depths on the level populations assumed to be at LTE, following the formula 


W, 




In(DhP) = In(^) - 




up 




9up 


Q 


rot 


kT ~ ^ 

A/j.rot 


( 3 ) 


where Wot is the total column density of the species in question, Wp is the observed column 
density of the upper state of the species with an upper energy including the opacity 
effect, Qrot is the partition function, is the beam solid angle, is the source solid angle. 
Cr is given by 


a = 


r 


1 — exp(—r) 

with r being the optical depth, r can be expressed as 


T = 


C ^ul fi^up.^tot 

AD Qrot 


exp( 


Eun \ / / hr'o 

)(exp(- 


kTr 


rot 


kE 


)-l) 


rot 


( 4 ) 

( 5 ) 


where c is the speed of light, Wi is the Einstein-A coefficient of spontaneous emission, AV 
is the FWHM hxed to 3 km/s. We performed a reduced Xred minimization by running a 
grid of 125,000 models covering a large parameter space in rotational temperature Trot (50 
values between 10 and 500 K), total column density in the source Wot (50 values between 
10^® and 10^° cm“^), and source size Og (50 values between 0.04 and 2 "). At LTE, the 
column density of every upper state Wp can be derived for each set of Wot, Wt, and source 
solid angle = 9"^ according to equation ([3]). The best-ht model populations are plotted 
together with the observed populations of the levels in Figures 191 and [TUI and are marked by 
red cross symbols. Tables [5] and [6] summarise the parameters of the best-ht models. 
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We started the PD analysis by simultaneously modelling the population distribution 
of ^^CH^OH and ^^CH.^OH. F or thi s purpose, we assumed a abundance ratio 

( 120021 ) and the same rotational temperature for the two 


of 70 following 


Boogert et al 


isotopologues. The population modelling of high-energy ^^CHsOH (400 < Eup < 1100 K) 
and low-energy ^^CHsOH optically thin transitions allowed us to constrain the rotational 
temperature through the determination of the slope of the level populations, leaving only 
a degeneracy between Wot and Og. Since the optical depth r of each level depends on the 
total column density Wot, low-energy optically thick transitions from ^^CHaOH can be used 
to constrain Wot and Og- 

Table [5] presents the results of the PD analysis of the methanol population distribution. 
The methanol emission is relatively well modelled by the PD model for the two sources. 
The PD analysis converges towards one single set of input physical parameters (Wot, 

^rot) with Xred ^f about 1.5-2 and with uncertainties up to 50 % at a one sigma level for the 
column densities. The PD model was able to reproduce the population of most transitions 
within their uncertainties except the population of some low-energy transitions of ^^CHsOH 
in IRAS2A and of ^^CHsOH in IRAS4A which tend to be underestimated. The rotational 
temperatures of methanol deduced from the PD analysis are similar in IRAS2A and 
IRAS4A (~ 140 K). However, the source size derived for IRAS2A is larger than for IRAS4A 
(~ 0.36" versus ~ 0. 20"). The source si ze deduced for IRAS2A is in good agreement with 


the size estimated by 


Maret et al 


(120041) but we found a smaller source size for IRAS4A, by 


a factor of 2.5. The source size of IRAS2A of 0.36 ", corresponding to a radius of 42 AU at 
23 5 pc, is also c o nsiste nt with the FWHM emission size of the CH 3 OH transitions deduced 


by 


Maury et al. 


(120141) from elliptical gaussian fits that range between 25 AU and 70 AU. 


For all other COMs, the low number of observed transitions and the relative high 
uncertainty on the derived column density of each level generates a degeneracy between 
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Table 5: Results of the PD analysis of the methanol emission. 

IRAS2A IRAS4A 


X?ed 

1.4 

2.1 

A^(CH 30 H) (cm-2) 

X 10^® 

1.6l°;® X 10^0 

NC^CHaOH) (cm-2) 

X 10^® 

2.31J1 X 10^7 

0s (") 


0.201^4 

Trot (K) 

140120 

1401” 


the input parameters, and prevents the PD model to converge towards one single set of 
input parameters: the observations are overhtted and can be reproduced by a large range 
of parameters giving lower than 1. Therefore, we decided to £x the source size for 
the COMs emission to the size of the methanol emission, assuming that all COMs will 
evaporate with methanol in spite of their slightly different binding energies. Even by hxing 
the source size, the analysis of the glycolaldehyde and ethyl cyanide populations towards 
the two sources and of the ethanol population towards IRAS4A did not allow us to converge 
towards one set of A^tot, and Trot- For other molecules, we were able to deduce a unique 
column density with relatively small uncertainties. As seen in Fig. |9]and[T0l the observed 
population distribution is well reproduced by the best £t model since most of the column 
densities of upper energy levels predicted by the best £t model he within the range of 
uncertainties of the observed data. The best £t models of the PDs generally consist in lower 
rotational temperatures and higher total column densities than the values derived with the 
RDs in order to reproduce the population of the optically thick transitions. For instance, 
the spread of the population distribution of the low upper energy levels (E^p < 120 K) 
of CH 3 OH, ^^CHsOH, HCOOCH 3 , or CH 3 CN are explained by opacity effects. For these 
species, transitions showing a decrease of their population are optically thick with r higher 
than 1 . 
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For species where the PD analysis was not able to converge towards one set of inpnt 
parameters (namely glycolaldehyde, ethanol, and ethyl cyanide) and for molecules showing 
only one detected transition, we hxed the source size and the rotational temperature to the 
values found for methanol. Results of this analysis are shown in Table [ 6 l 


Most transitions of COMs whose collision rates have been computed (CH 3 OH, CH 3 CN, 
and HC 3 N) have critical densities that range between 10^ and 10^ cm“^ at 100 K. They are 
therefore likely lower than the densities found in the hot corinos of IRAS2A and IRAS4A 
{nm > 1-3 X 10^ cm“^ following the density profiles of the two envelopes derived by 


Maret et af 


20041) . Given the good fit to the observational data with our LTE PD analysis. 


it is likely that the observed species are at LTE. Most of the scattering of the population 
distribution can therefore be attributed to opacity effects only. 


5. Discussion 

5.1. Abundances in IRAS2A and IRAS4A 

Methanol is likely the most abundant complex organics, and is believed to be the 
precursor molecule of several COMs. It is therefore worth comparing the abundance of 
COMs with respect to methanol to quantify the efficiency of their formation. Moreover, 
column densities of COMs and methanol have been derived with similar methods and from 
the same observational data, the estimates of the abundance ratios are therefore more 
accurate than the absolute abundances derived with respect to H 2 . Tables H] and [H] list 
the abundance ratios of the COMs with respect to methanol for the RD and PD analyses. 
The two targeted sources seem to have a similar chemical composition since the COMs 
abundance ratios differ by only a factor of two at maximum. Methyl formate is the most 
abundant COM of our sample, with an abundance of 1.5-3 % followed by ethanol (1-1.5 %) 
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Table 6: Results from the population diagram analysis for the COMs emission. 


Molecule 

Nhc 

Trot 

Source size 

A“ 

^meth 


(cm“^) 

(K) 

(") 



IRAS2A 

CH 3 OH 

S.OtiJ X IQi® 

140t^2o 

U.OO_o.o4 

l.Otn X 10 -® 


13 CH 3 OH 

7 . 1 ^ 2 ;^ X IQi® 

140t^2o 

U.OO_o.o4 

1.4t°;® X 10"® 


HCOOCH 3 

X IQi® 

160t®o 

0.36 ^ 

1.6t°;® X 10-® 

1 6 +^-^ 

CH 3 CN 

2.0^11 X IQi® 

130t^®° 

0.36 ^ 

4.0to;8 X 10 ”® 

^•^^-0.25 

CH 3 OCH 3 

X IQi® 

not®® 

0.36 ^ 

l.Ot!^;^ X 10 -® 

1 O"'"®!' 

C 2 H 5 OH 

X 10 ^® 

r)'yi^+230 

^<U_8o 

0.36 

1 . 6 t!^J X 10 "® 

1 6 +^-^ 
1 . 0_1 2 

HCOCH 2 OH 

6.8 X 10^® 

140^^ 

0.36 

1.4 X 10"® 

0 14+0.05 

U.l^_0.08 

C 2 H 5 CN 

1.5 X 10^® 

140^^ 

0.36 

3.0 X 10-1® 

0 

0 

CO 

0 

1 + 

0 0 

b b 

HC 3 N 

9.3 X lOi-* 

IdO'^ 

0.36 

1.9 X 10-1® 

0 01 

H 2 I 3 CO 

6.3 X IQi® 

140^^ 

0.36 

1.3 X 10“® 

0 1 Q+0.05 

U.lOoo? 

NH 2 CHO 

1.2 X IQi® 

140^^ 

0.36 ^ 

2.4 X 10“® 

0 

^•^^-0.14 

CH 2 CO 

6.8 X 10^® 

140^^ 

0.36 ^ 

1.4 X 10“® 

0 14+005 

U.i^0.08 

IRAS4A 

CH 3 OH 

1.6t°;® X IQi® 

140t®o 

0 

U-^^_0.04 

4.3t^;® X 10-1' 


13 CH 3 OH 

2.3ti;^ X 10 ^^ 

140t®o 

U.ZU_o.o 4 

6.2t®;® X 10-® 


HCOOCH 3 

5 . 0 t^;^ X IQi^ 

6 o;?o 

0.20 ^ 

1.4tJ;^ X 10-® 

Q 1 +3.5 

0 .i_ 2 .i 

CH 3 CN 

6 . 3 t?;^ X IQi® 

200t“° 

0.20 ^ 

1.7tJ;^ X 10-® 

0 QQ+0.30 
U.OJ_o ,24 

CH 3 OCH 3 

1 . 6 t!^;^ X IQi^ 

801^0 

0.20 

4.3tg;| X 10-® 

1 n'*'®® 

C 2 H 5 OH 

1.6 X 10^^ 

140^^ 

0.20 

4.3 X 10-® 

1-0+:^ 

HCOCH 2 OH 

4.8 X 10^® 

IdO'^ 

0.20 

1.3 X 10-® 

O.SOt®;!® 

C 2 H 5 CN 

6.4 X IQi® 

140^^ 

0.20 ^ 

1.7 X 10-1® 

0.040t°;°2o 

HC 3 N 

2.9 X 10^® 

140^^ 

0.20 ^ 

7.8 X 10-11 

0-018t°:°?o 

H2I3CO 

1.1 X IQi® 

140^^ 

0.20 ^ 

3.0 X 10-1® 

u.uoyo.o 40 

NH 2 CHO 

1.9 X 10^® 

IdO'^ 

0.20 

5.1 X 10-1® 

0 

^■-‘-^- 0.07 

CH 2 CO 

3.4 X 10^® 

IdO'^ 

0.20 

9.2 X 10-1® 

0 2 i+'^ io 

u.zi_o 12 


Note. — a: The abundances relative to H2 are obtained from fV(H2) derived at 145 GHz in Table [T] 
assuming an homogeneous H2 column density within the beam, b: The abundances relative to CH3OH 
were computed from A^(CH30H) derived from the PD analysis of the CH3OH and ^^CH30H emissions and 
adapted for the same source size, c: The source size was assumed to be equal to that of methanol when the 
size could not be constrained, d: The source size and rotational temoeratures were assumed to be eoual to 
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and di-methyl ether (1 %). Other COMs are detected with abundances lower than 1 %: 
glycolaldehyde and methyl cyanide show abundances of 0.15 and 0.40 % respectively whilst 
ethyl cyanide is detected with an abundance of 0.03-0.04 %. 


Table 0] compares the column densities and abundance r atios c 


analysis with previous single-dish studies by 


Bottinelli et ah 


Bottinelli et ah 


il2004ah : 


educed from the RD 


Maret et ahl (2005); 


(120071) carried out towards 1RAS2A and 1RAS4A. The column densities 


obtained from these previous observations suffer from several limitations: most detected 
transitions have low upper energy levels and the large beam of single-dish telescopes 
encompasses the cold envelope where weakly excited lines may have contaminated the 
hot corino emission (see the interferometric maps in Figures |6] and [7]). Consequently, 
the rotational temperatures and column densities derived with the PdBi are higher (Tj-ot 
= 80-290 K and 300-360 K in this work towards 1RAS2A and 1RAS4A versus 100 and 
25 K respectively and higher column densities up to one order of magnitude), since the 
interferometric observations probe material closer to the central protostars. The abundance 
ratios deduced from previous single-dish studies have also higher relative uncertainties than 
in this work due to the different telescope calibrations since the observations of methanol 
and COMs have been carried out separately. The abundances relative to methanol derived 
from our interferometric observations therefore differ from the abundances obtained with 
the single-dish observations, the latter being usually overestimated. The methyl formate 
abundance deri v ed in 1RAS2A of 2 % is consistent with the upper limit of 85 % by 

(120071 ). However , the abundance derive d in 1RAS4A of ~ 3 % is 18 times 


Bottinelli et al. 


lower than the valu e derived by 
CH3OH derived by Maret et al. 


Bottinelli et al. 




(j2004al ) and using the column density of 


(120051) . The higher abundance of methyl formate with 


respect to methanol in 1RAS4A derived from single-dish o 


lower column density of methanol derived in 


Maret et al 


Dservations is explained by the 


( 20051 ) assuming optically thin 


emission. For the same reasons, the abundances of methyl cyanide derived in the two 
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sources by our PD analysis are a. 


obtained by 


Bottinelli et al.l ( 2007 


so lower, by a factor of 3 to 6, than the abundances 


). Di-methyl ether, ethanol and ethyl cyanide have not 


been detected with single-dish telescopes towards IRAS2A and IRAS4A but their upper 
limits agree well with our observations. 

We report here the first detection of glycolaldehyde in low-mass protostars other 
than IRAS 16293. Glycolaldehyde co-exists with its isomer methyl formate with a 
[HC00CH3]/[HC0CH20H] abundance ratio of 12^2 towards IRAS2A and of 10l4° towards 


IRAS4A. These abundance ratios are similar to the ratios of ~ 13 founc 
and B of the IRAS 16293 protostellar binary system by 


Jprgensen et al. 


in the sources A 


(120121) from high 


angular resolution ALMA observations. They are a 


towards SgrB2(' N') ranging from 52 in the hot core (IHollis et al.l 120011) to 5 found on more 


extended scales (IHollis et al. 


2000 |). 


so consistent with the ratios derived 


5.2. Comparison with Other Sources 


The COMs abundances derived in this work are compared with other published data 
of low-mass, intermediate-mass, and high-mass hot cores obtained with single dish and 
interferometric sub-mm telescopes in Figure [8] showing the COMs abundances as a function 
of the protostar luminosity. Table [T7] lists the abundances of selected COMs towards the 
hot cores shown in Fig. [H] along with the references. For this purpose, we only selected 
observational studies where several transitions of methanol were detected. For most of the 
works, COMs and methanol were detected simultaneously and the abundance ratios were 


and LTE population ( 

MacDonald e 

al. 

1996 

Ikeda et al. 2001: 

Beuther et al. 2007, 2009 

Palau et al. 

2011 


Oberg et al. 

2011 

, 

201 

4), t 

V neglecting optic 

ally thick lines in the RD 

g the opacity of the lines 

analysis ( 

Bisschop et al. 

2007; 

Isokoski et al. 

2013) or bv takin 
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into account in t heir model (|Niimmelin et al 


Neill et al 


201J). 


Fuente et al 


2000 


Qin et al. 


2010 


Crockett et al. 


2014 : 


(120141) derived the methanol abundances from the ^^CHaOH 


isotopologue. 


Schoier et ah 


or I 


4.AS 16293, we combined the CH 3 OH abso 


(120021) with the COMs abundances obtained by 


ute abunda n ce de rived by 


Jaber et ah 


(120141) ■ both 


from a radiative transfer model 
glycolaldehyde abundance from 


ing, to obtain t 


le abundance ratios. We also derived the 


Jorgensen et ah 


( 120121 ) . 


Overall, the abundance ratios of COMs estimated towards IRAS2A and IRAS4A in 
this work tend to be lower than the abundances derived in other low-mass protostars from 
single-dish observations. For example, methyl formate and methyl cyanide have been 
detected towards 5 other low-mass sources and show abundances of 5-50% and 0.8-1.7% 
respectively, representing a factor of 1.5-15 and 1.5-3 higher than in our work. Since 
these observations suffer from the same limitations than the single-dish observations by 


Bottinelli et al. 


(120043, 


20071) presented in the previous section (low number of detected 


transitions, low upper energy levels of detected transitions, large beam encompassing the 
external envelope), the discrepancy likely comes from the differences in the observational 
methods and does not necessarily reflect differences in the chemistry between the sources. 


Table 171 summarises the mean abundance ratios of COMs in low-mass, intermediate- 
mass and high-mass protostars. Along with Fig. IHl they allow us to investigate any possible 
correlation of the COMs abundances with the protostar luminosity. For each molecule, 
the data has been htted by a linear curve depicted by the dashed curve whose slope a is 
shown at the top left of each panel of Fig. IHl by considering detected ratios only. It can be 
noticed that the abundances of the six COMs tend to slightly increase with the protostar 
luminosity. However, for all species but C2H5CN, the increase remains negligible compared 
to the dispersion of the abundance ratio values. We can conclude that the abundance ratio 
of these COMs stays relatively constant with the protostar luminosity within six orders of 
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magnitude. In spite of their lower luminosities, inducing lower temperatures, and smaller 
sizes, low-mass protostars seem to be as chemically complex as high-mass protostars. 

Table 7: Averaged COMs abundances (in %) with respect to methanol. 


Molecule 

LMP 

IMP 

HMP 

G13 

ROl 

HCOOCH 3 

141^' 

1.3 ±0.2 


0.082 

- 0.84 

0.013 ■ 

- 0.67 

CH 3 OCH 3 


1.2 

sotll 

0.44 - 

- 0.74 

0.62 

- 18 

HCOCH 2 OH 

0.381°:?^ 

/ 

3.1 

0.54 

- 1.3 

- 


C 2 H 5 OH 

1 q+0.3 
-‘-•'^-0.3 

q 9-r5.9 

'J-^-2.0 

6.4+®;? 

0.50 

- 2.6 

15 - 

23 

CH 3 CN 

-‘-•■‘--0.5 

0.52 


0.034 

- 0.45 

3.0 - 

3.7 

C 2 H 5 CN 

u.uau_o.o55 

/ 


0.052 

- 0.79 

- 



Note. — LMP, IMP, and HMP stand for Low-Mass (L < lOOL©), Intermediate-Mass (lOOL© < L < 
IO^Lq), and Higli-JVIass ^ L) Protostars. Tire two abundances in the column G13 are the peak 


abundances of the “Fast” and “Slow” models of 


GarrodI ( 20131) . The two abundances in the column “ROl” are 


the pe ak abundances of the models including ammonia at T = 100 and T = 300 K of 

(I2OO1I). 


Rodgers and Gharnlev 


5.3. Comparison with Chemical Model Predictions 


We compare the observed abundance ratios in low-mass, intermediate-mass, and 
high-mass pro tostars with the results of two astrochemical models in Table [3 The model of 
GarrodI (120131) is a multilayer gas-grain astrochemical model in which COMs are assumed to 
be mostly formed at the surface of interstellar ices during the warm-up protostellar phase 
(30 K < T < 100 K) through radical recomb ination induced by the UV p hotodissociation 


of the main ice components. In the model of 


Rodgers and CharnlevI fl 200 ll) . COMs are only 


formed by warm gas phase chemistry for a set of constant physical parameters representative 
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of hot cores (uh = 10^ cm“^; T = 100 and 300 K) after the sublimation of interstellar ices 
with typical ice composition, already containing C 2 H 5 OH with an abundance of 20 % with 
respect to CH 3 OH. Although the rates of some key reactions for the formation of COMs 
have been lowered meanwhile, such as the methyl cation transfer reaction between H 2 CO 
and CHaOH^ or the electronic recombination of protonated COMs ions, this model still 
provides a good basis to estimate the formation efficiency of methyl formate, di-methyl 
ether and methyl cyanide in the gas phase. 


Both models tend to underpredict the abundance of methyl formate relative to 
methanol observed in lo w-mass and hig h-mass protostars by at least 1-2 orders of magnitude. 


Moreover, the model of 


GarrodI ( 2013 1 also underpredicts the [HC 00 CH 3 ]/[HC 0 CH 20 H] 


abundance ratio by two orders of magnitude since it seems to reproduce well the observed 
abundance of glycolaldehyde. This comparison suggests that the chemical network forming 
methyl formate either in the gas phase or on ices is still incomplete. Possible alternative 
branching ratios for the photodissociation of CH 3 OH on ices or other gas phase reactions 
involving HCOOH could enhance its formation. The observed abundances of di-methyl 
ether and methyl cyanide with r e spect to methanol can be reproduced by the gas phase 


model of iRodeers and CharnlevI fl200lf l only. The absolute abundances of these two 
molecules are similar in the two models ( 10 “® — 10 “^ for di-methyl ether and 10 “® — 10 “® for 
methyl cyanide), showing that warm gas phase chemistry tends to be as efficient as surface 
chemistry to produce these COMs. However, the difference in the a bundances comes from 
the effi cient destruction of methanol in the warm gas in the model of 


Rodgers and Charnley 


(120011 ) increasing the abundance ratio of COMs. No efficient formation routes in the 


gas phase have been proposed for glycolaldehyde and ethanol but their formation at the 
surface of interstellar ices seem to be efficient enough to reproduce the observations towards 
low-mass and high-mass protostars. The abundance of ethyl cyanide shows an increase 
of almost two orders of magnitude between low-mass and high-mass protostars. Grain 
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surface chemistry is able to reproduce the observations towards low-mass protostars but not 
towards high-mass hot cores. It is also possible that models also missed gas phase reactions, 


as it was the case for methy 


recognised by 


Balucani et al 


brma te, where a new gas phase reaction has been recently 


fcnibh . 


6. Conclusions 

In this work, we have presented interferometric multi-line observations of methanol 
(CH 3 OH, ^^CHsOII) and various COMs towards the two bright low-mass protostars 
NGC1333-IRAS2A and -IRAS4A, carried out with the Plateau de Bure interferometer at 
an angular resolution of ~ 2". We summarize here the main conclusions of this work: 

- Except for methanol, none of the transitions from complex organics are spatially resolved 
with the interferometer, conhrming that most of the emission comes from the inner arcsec 
from the central protostars. 

- Low-energy transitions {E^p < 120 K) of methanol show extended emission around 
IRAS4A and trace the outflows driven by the central protostar. 

- We detected a high number of transitions (up to 45 for methanol) from COMs with a wide 
range of upper energy levels (up to 1000 K for methanol) allowing us to accurately derive 
their column densities with the use of LTE methods. 

- The abundances derived in this work, of a few percents for methyl formate and di-methyl 
ether and of ~ 0.5% for methyl cyanide for instance, seem to be slightly lower than the 
abundances found towards other low-mass protostars. However, the difference likely comes 
from differences in the observational methods rather than different chemistries taking place 
in these protostars. 

- We report the hrst detection of glycolaldehyde in low-mass protostars other than IRAS 
16293, with a methyl formate-to-glycol aldehyde abundance ratio of 12 and 10 in IRAS2A 
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and IRAS4A, respectively. 

- The analysis of the variation of the COMs abnndance ratios with the protostellar 
Inminosity shows that low-mass hot corinos seem to be as chemically rich as high-mass hot 
cores, in spite of their lower temperatnres and their smaller sizes. 

- Comparison with theoretical models shows that the two theories of COMs formation, 
either in the gas phase or at the snrface of interstellar ices, tend to nnderpredict the 
formation of methyl formate by one to orders of magnitnde. The comparison of the 
abnndance ratios of other molecnles favonrs a gas phase formation for di-methyl ether and 
methyl cyanide and a formation on ices for ethanol, ethyl cyanide, and glycolaldehyde. 

The anthors are gratefnl to the anonymous referee whose comments contributed to 
improve the quality of the present paper. This work was supported by NASA’s Origins 
of Solar Systems and Exobiology Programs. V. T. acknowledges support from the NASA 
postdoctoral program. A. L-S and C. C acknowledge hnancing from the french space agency 
CNES. 
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Fig. 5.— Top) Deconvolved FWHM linewidths of CH 3 OH, i^CHgOH, HCOOCH 3 , and 
CH 3 CN, obtained from gaussian fits of the WideX spectra at the source positions, as a 


function of E^p. Bottom) Comparison of the FWHM linewidths deduced from the high 
resolution spectra (plus signs with solid error bars) with the deconvolved FWHM linewidths 
from the WideX spectra (cross signs with dotted error bars) for the lines observed with the 
high spectral resolution correlators. 










































































































35 



RA offset (arcsec) 


Fig. 6 .— Integrated maps of selected lines of complex organic molecules towards IRAS2A. 
For each species where several transitions have been detected, a low-excited and a high- 
excited transition are shown. Contour levels are in step of Sa. la rms noise levels are the 
following: 38.5 and 18.2 mJy km/s for CH 3 OH, 13.3 and 28.9 mJy km/s for ^^CHsOH, 13.0 
and 14.6 mJy km/s for for IICOOCH 3 , 35.3 and 26.5 mJy km/s for CII 3 CN, 15.5 and 9.64 
mJy km/s for CII 3 OCH 3 , 11.6 and 14.5 mJy km/s for C 2 II 5 OII, 21.7 and 19.5 mJy km/s 
for HCOCH 2 OH, 10.9 and 10.8 mJy km/s for C 2 H 5 CN, 24.8 mJy km/s for IIC 3 N, 13.6 mJy 
km/s for H^^CO, 17.6 mJy km/s for NII 2 CHO, 12.0 mJy km/s for CII 2 CO. The direction of 
the outflow is depicted by the dotted line. 
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RA offset (arcsec) 


Fig. 7.— Integrated maps of selected lines of complex organic molecules towards IRAS4A. 
For each species where several transitions have been detected, a low-excited and a high- 
excited transition are shown. Contour levels are in step of Sa. la rms noise levels are the 
following: 81.9 and 22.8 mJy km/s for CH 3 OH, 27.4 and 32.9 mJy km/s for ^^CHsOH, 25.4 
and 24.5 mJy km/s for for IICOOCH 3 , 43.6 and 24.0 mJy km/s for CII 3 CN, 16.9 and 11.4 
mJy km/s for CII 3 OCH 3 , 12.4 and 19.8 mJy km/s for C 2 II 5 OII, 18.8 and 17.2 mJy km/s 
for HCOCH 2 OH, 11.6 and 16.0 mJy km/s for C 2 H 5 CN, 109 mJy km/s for IIC 3 N, 27.8 mJy 
km/s for H^^CO, 14.3 mJy km/s for NII 2 CHO, 22.4 mJy km/s for CII 2 CO. The direction of 
the outflow is depicted by the dotted line. 
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Fig. 8.— Abundances of COMs detected in several transitions with respect to CH 3 OH 
towards a sample of low-mass (L < IOOLq; blue), intermediate-mass (IOOLq < L < IO^Lq; 
green), and high-mass (L > IO^Lq; red) hot cores as function of the protostar luminosity. 
Cross and plus symbols represent abundance ratios derived with single-dish telescopes and 
interferometers, respectively. The abundance ratios derived in this work with the Population 
Diagram analyses are denoted by squares with error bars. The sample of hot cores is listed 
in Table [IT] along with their physical and chemical properties. 
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Table 8 . Line parameters of CH 3 OH lines observed towards IRAS2A and IRAS4A-NW. 


IRAS2A IRAS4A 


N 

Frequency 

Transition 

Eitp 

^ul 

Beam size 

Source size “ 

Flux ^ 

dVw "" 

Beam size 

Source size (PA) 

Flux ^ 

dVw ^ 


(GHz) 


(K) 


("x".-’) 

("x". °) 

(Jy km/s) 

(km/s) 

("x". p 

("x". p 

(Jy km/s) 

(km/s) 

1 

145.093754 

3o-2o El, ut—0 

27.1 

1.23(-5) 

2.08x1.65 (30) 

1.65x1.04 (19) 

2.62±0.53 

- 

2.16x1.73 (25) 

outflow 

2.54±0.51 

- 


145.097435 

3o-2q E2, ut—0 

19.5 

1.10(-5) 










145.103185 

3o-2o A”*", ut—0 

13.9 

1.23(-5) 









2 

165.050175 

Ii-Iq El, ut—0 

23.4 

2.35(-5) 

2.31x1.72 (110) 

1.01x0.87 (11) 

1.28±0.26 

5.8 

2.39x1.77 (114) 

outflow 

0.38±0.08 

6.2 

3 

165.061130 

2i-2q El, ut=0 

28.0 

2.34(-5) 

2.31x1.72 (110) 

1.18x0.83 (50) 

1.49±0.30 

5.2 

2.39x1.77 (114) 

outflow 

0.56±0.12 

7.0 

4 

143.865795 

3i-2i A"*”, ut=0 

28.3 

1.07(-5) 

2.24x1.77 (25) 

0.83x0.63 (10) 

1.02±0.21 

6.7 

2.24x1.77 (25) 

outflow 

0.44±0.13 

7.5 

5 

165.099240 

3i-3o El, ut=0 

35.0 

2.33(-5) 

2.31x1.72 (110) 

1.19x1.07 (30) 

1.44±0.29 

6.0 

2.39x1.77 (114) 

outflow 

0.39±0.08 

5.2 

6 

145.124332 

3o-2o A“ , ut—0 

51.6 

6.89{-6) 

2.08x1.65 (30) 

1.05x0.71 (32) 

3.18±0.64 

- 

2.16x1.73 (25) 

outflow 

1.65±0.33 

- 


145.126191 

32-22 El, ut—0 

36.2 

6.77(-6) 










145.126386 

32-22 E2, ut—0 

39.8 

6.86(-6) 










145.131864 

3i-2i El, ut=0 

35.0 

1.12(-5) 










145.133415 

32-22 A”*", vt—0 

51.6 

6.89(-6) 









7 

165.190475 

4i-4q El, ut=0 

44.3 

2.32(-5) 

2.31x1.72 (110) 

1.26x1.02 (76) 

1.96±0.40 

6.0 

2.39x1.77 (114) 

outflow 

0.64±0.13 

9.5 

8 

165.369341 

5i-5o El, vt^O 

55.9 

2.31(-5) 

2.54x1.71 (113) 

1.10x0.92 (42) 

1.59±0.32 

6.1 

2.38x1.76 (114) 

outflow 

0.68±0.14 

5.5 

9 

165.678649 

61-60 El, ut=0 

69.8 

2.30(-5) 

2.54x1.71 (113) 

1.08x0.99 (54) 

1.52±0.31 

5.9 

2.38x1.76 (114) 

outflow 

0.551^0.11 

0.543tA-ll 

6.3 

10 

166.169098 

7i-7o El, ut—0 

86.1 

2.28{-5) 

2.54x1.71 (113) 

1.06x0.99 (30) 

1.46±0.29 

5.8 

2.38x1.76 (114) 

outflow 

5.8 

11 

166.898566 

81-80 El, ut=0 

104.6 

2.28(-5) 

2.54x1.71 (113) 

1.02x0.56 (44) 

1.45±0.29 

5.7 

2.38x1.76 (114) 

outflow 

0.51zd0.11 

5.0 

12 

143.169517 

73-82 El, ut—0 

112.7 

4.13(-6) 

2.31x1.82 (26) 

point 

0.63±0.22 

6.1 

2.31x1.82 (26) 

outflow 

0.24±0.08 

5.8 

13 

144.728359 

32-22 A+. vt = 1 

312.6 

6.78(-6) 

2.08x1.65 (30) 

0.68x0.64 (-73) 

0.80±0.16 

6.9 

2.16x1.73 (25) 

point 

0.23±0.05 

- 


144.728782 

32-22 E2, vt — 1 

378.5 

6.83(-6) 










144.729074 

32-22 A“ , ut —1 

312.6 

6.78(-6) 









14 

144.733262 

32-22 El, vt — 1 

413.8 

6.80(-6) 

2.08x1.65 (30) 

0.38 

0.96±0.19 

9.4 

2.16x1.73 (25) 

0.48 

0.42±0.09 

7.6 


144.734429 

3i-2i El, ut = l 

305.4 

1.09(-5) 










144.736349 

3o-2o El, Vt — 1 

314.5 

1.22(-5) 









15 

144.589854 

3i-2i A+, vt^l 

339.1 

1.09(-5) 

2.08x1.65 (30) 

0.59 

0.37±0.08 

6.0 

2.16x1.73 (25) 

0.3 

0.21±0.06 

6.7 

16 

144.878576 

3i-2i A“ , ut —1 

339.2 

1.09(-5) 

2.08x1.65 (30) 

0.35 

0.43±0.09 

6.7 

2.16x1.73 (25) 

0.91x0.75 (-5) 

0.23±0.06 

5.6 

17 

143.108385 

17o-17_i E, vt=0 

366.8 

6.44(-6) 

2.06x1.67 (26) 

0.61x0.50 (0) 

0.52±0.11 

6.9 

2.06x1.67 (26) 

0.50 

0.18±0.04 

6.9 

18 

166.773281 

5i-5o A"*", vt = l 

390.0 

1.85(-5) 

2.54x1.71 (113) 

0.91x0.46 (67) 

0.90±0.19 

6.4 

2.38x1.76 (114) 

0.58 

0.45±0.10 

7.5 

19 

165.074355 

146-155 El, ut-0 

422.4 

5.36(-6) 

2.31x1.72 (110) 

0.90x0.62 (37) 

0.47±0.10 

5.5 

2.39x1.77 (114) 

2.02x1.01 (-53) 

0.30±0.06 

6.0 

20 

144.750264 

3i-2i E2, vt = l 

427.3 

1.08(-5) 

2.08x1.65 (30) 

point 

0.28±0.06 

7.2 

2.16x1.73 (25) 

0.15 

0.12±0.04 

6.7 

21 

144.768193 

3o-2o A+, vf — l 

437.5 

1.22(-5) 

2.08x1.65 (30) 

0.46x0.40 (49) 

0.28±0.06 

6.8 

2.16x1.73 (25) 

1.23 

0.09±0.03 

4.4 

22 

144.572025 

3o-2o A”*", vt—2 

522.1 

1.42(-5) 

2.08x1.65 (30) 

0.52 

0.21±0.04 

6.3 

2.16x1.73 (25) 

1.61x1.03 (-60) 

0.15±0.04 

5.5 


144.571262 

32-22 E2, vt—2 

658.8 

6.74{-6) 









23 

144.583961 

3i-2i E2. i>t=2 

545.9 

1.09(-5) 

2.08x1.65 (30) 


0.20±0.04 

4.0 

2.16x1.73 (25) 


0.15±0.05 

16.3 

24 

166.843662 

II2-IO3 El, vt = l 

553.0 

1.19(-6) 

2.54x1.71 (113) 

0.48 

0.07±0.02 

3.5 

2.38x1.76 (114) 

- 

< 0.04 

- 



Table 8—Continued 


25 

164.299104 

152-14i A”, vt = l 

576.0 

1.44(-5) 

2.31x1.72 (110) 

0.94x0.64 (-86) 

0.43±0.10 

5.5 

2.39x1.77 (114) 

1.95x0.94 (-83) 

0.21±0.05 

3.6 

26 

142.276432 

97 -IO 6 E2, vt = l 

627.5 

l.ll(-6) 

2.06x1.67 (26) 

1.66x0.50 (-42) 

0.03±0.01 

6.6 

2.06x1.67 (26) 

0.41 

0.02±0.01 

4.8 

27 

144.281736 

3i-2i A"^, v±=2 

696.6 

1.07(-5) 

2.08x1.65 (30) 

point 

0.08±0.02 

7.2 

2.16x1.73 (25) 

0.78 

0.06±0.03 

5.3 

28 

144.530553 

3o-2o El, ut—2 

748.1 

1.22(-5) 

2.08x1.65 (30) 

0.34 

0.04±0.01 

4.0 

2.16x1.73 (25) 

1.15 

0.05±0.01 

4.1 

29 

144.499723 

3i-2i El, v±=2 

812.8 

1.05(-5) 

2.08x1.65 (30) 

0.96 

0.03±0.01 

7.0 

2.16x1.73 (25) 

point 

0.03±0.02 

5.7 

30 

164.486238 

262"26i E2, vt=0 

843.0 

2.69(-5) 

2.31x1.72 (110) 

0.85x0.43 (38) 

0.30±0.07 

5.8 

2.39x1.77 (114) 

1.25 

0.26±0.05 

5.8 

31 

144.195014 

246-255 E2, vt=0 

884.4 

5.16(-6) 

2.08x1.65 (30) 

1.76x0.69 (-59) 

0.05±0.01 

6.0 

2.16x1.73 (25) 

point 

0.07±0.03 

5.9 

32 

144.437702 

lOg-lOaO E2, -ut —1 

922.9 

2.82{-6) 

2.08x1.65 (30) 


< 0.05 


2.16x1.73 (25) 

1.39 

0.04±0.01 

7.6 

33 

144.807264 

119 -llaO El, ^t-1 

948.3 

4.94(-6) 

2.08x1.65 (30) 


< 0.05 


2.16x1.73 (25) 

1.39x0.67 (-2) 

0.07±0.02 ' 
0.10±0.04C^ 

6.9 

34 

145.068727 

195-206 E2, vt = l 

985.9 

8.88{-6) 

2.08x1.65 (30) 

1.33 

0.07±0.02 

6.5 

2.16x1.73 (25) 

0.6 

7.0 

35 

163.526070 

284-275 A“, vt=0 

1021.9 

8.13(-6) 

2.31x1.72 (110) 

1 

0.04±0.02 

4.7 

2.39x1.77 (114) 

point 

0.10±0.04 1 

7.0 

36 

164.554640 

284-275 A+, vt=0 

1021.9 

8.28(-6) 

2.31x1.72 (110) 

2.57x1.42 (28) 

0.08±0.03 


2.39x1.77 (114) 

point 

0.04±0.02 

5.7 


Note. — a: Size of the best gaussian fit to the visibilities. Gaussian fits were performed on the channel showing the peak emission. When an 
elliptical gaussian fit was not successful, a circular gaussian was attempted. Gaussian fits resulting in a size of 0 ^^are marked by “point”. Lines 
dominated by a molecular outflow are marked by “outflow”. See text for more details. 

b: Flux derived from a circular mask with a diameter equal to the major axis of the beam of the telescope given for each transition in the Table. 
The errors on the fluxes were computed as the quadratic sum of the statistical error and the calibration uncertainty (about ~ 20 %). 
c: Non-deconvolved FWHM linewidths derived from a gaussian fit of the WideX spectra towards the source positions. 




N 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 
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Table 9. Same as Table [8] but for ^^CHsOH. 


IRAS2A IRAS4A 


Frequency 

Transition 

Eup 

A„i 

Beam size 

Source size 

Flux 

dVw 

Beam size 

Source size 

Flux 

dVwr 

(GHz) 


(K) 

(s-b 

("x". °) 

("x". °) 

(Jy km/s) 

(km/s) 

C'x", 

("x". °) 

(Jy km/s) 

(km/s) 

141.603710 

3o-2o, A"^, ut—0 

13.6 

1.15(-5) 

2.06x1.67 (26) 

1.26x0.84 (-40) 

0.1424=0.032 

4.6 

2.06x1.67 (26) 

2.15x1.13 (-33) 

0.0994=0.026 

6.9 

165.575639 

2i-2q , El, vt =0 

27.9 

2.37(-5) 

2.54x1.71 (113) 

0.48 

0.1674=0.044 

5.2 

2.38x1.76 (114) 

2.27x1.39 (-62) 

0.1674=0.037 

3.6 

142.807657 

3i-2i A , vt=0 

28.3 

1.05(-5) 

2.06x1.67 (26) 

1.20x0.78 (18) 

0.1124=0.026 

6.6 

2.06x1.67 (26) 

point 

0.0694=0.020 

5.6 

165.609427 

3i-3o El, i>t=0 

34.6 

2.36(-5) 

2.54x1.71 (113) 

0.56 

0.2934=0.073 

6.4 

2.38x1.76 (114) 

0.56 

0.19i4o.042 

5.6 

165.690996 

4i-4o El, vt=0 

43.7 

2.34(-5) 

2.54x1.71 (113) 

1.04x0.71 (83) 

0.270=b0.068 

5.0 

2.38x1.76 (114) 

0.91 

0.172^.040 

5.0 

165.851224 

5i-5o El, vt^O 

55.0 

2.33(-5) 

2.54x1.71 (113) 

0.5 

0.294=b0.067 

5.0 

2.38x1.76 (114) 

2.43x0.73 (-55) 

0.1654^0.036 

4.1 

142.173740 

52-6i E2, vt=0 

60 

4.85(-6) 

2.06x1.67 (26) 

1.25 

0.0334=0.010 

6.4 

- 

- 

- 

- 

166.128782 

6i-6o El, vt=0 

68.6 

2.32(-5) 

2.54x1.71 (113) 

0.99x0.42 (25) 

0.3274=0.078 

5.2 

2.38x1.76 (114) 

1.74x0.75 (-88) 

0.1664=0.041 

5.5 

163.872900 

7q-6i El, vt=0 

76.5 

9.91(-6) 

2.31x1.72 (110) 

point 

0.156=b0.051 

5.4 

2.39x1.77 (114) 

1.80x0.74 (-45) 

0.1514=0.037 

4.5 

166.569486 

7i-7q El, vt—0 

84.5 

2.31(-5) 

2.54x1.71 (113) 

0.99x0.95 (-68) 

0.419=b0.093 

6.6 

2.38x1.76 (114) 

1.60x1.11 (-60) 

0.2334=0.050 

5.6 

142.896760 

62-7i A“ , ut—0 

85.4 

5.31(-6) 

2.06x1.67 (26) 

2.24 

0.0854=0.021 

8.1 

2.06x1.67 (26) 

1.45x0.77 (0) 

0.0904=0.023 

7.5 

165.280537 

62-7i A+, vt^O 

85.4 

7.77(-6) 

2.54x1.71 (113) 

point 

0.0674=0.020 

4.4 

2.38x1.76 (114) 

0.76 

0.0794=0.018 

8.4 

164.531587 

13i-122 A“, 'i;t=0 

222.3 

9.96(-6) 

2.31x1.72 (110) 

1.01 

0.1214=0.047 

3.6 

2.39x1.77 (114) 

2.15 

0.1054=0.028 

4.2 








Table 10. Same as Table [H] but for HCOOCH 3 . 


N 

Frequency 

(GHz) 

Transition 

Eup 

(K) 

A„, 

(s“b 


IRAS2A 




IRAS4A 



Beam size 

("x", “) 

Source size 

("x", °) 

Flux 

(Jy km/s) 

dV\Y 

(km/s) 

Beam size 

("x", °) 

Source size 

("x", °) 

Flux 

(Jy km/s) 

dV\^ 

(km/s 

1 

141.652995 

E, ll2,9"102,8 

43.2 

4.06(-5) 

2.06x1.67 (26) 

point 

0.090±0.024 

7.2 

2.06x1.67 (26) 

1.65x0.50 (-14) 

0.112±0.029 

7.0 

2 

141.667012 

A, ll2,9"102,8 

43.2 

4.06(-5) 

2.06x1.67 (26) 

0.2 

0.069±0.018 

7.2 

2.06x1.67 (26) 

0.93 

0.119±0.027 

6.6 

3 

143.234201 

E, 12 i,ii-lli,io 

47.3 

4.22(-5) 

2.31x1.82 (26) 

point 

0.073±0.028 

5.8 

2.31x1.82 (26) 

0.44 

0.091±0.035 

4.1 

4 

143.240505 

A, 12 i,ii-lli^io 

47.3 

4.23(-5) 

2.31x1.82 (26) 

point 

0.072±0.023 

5.8 

2.31x1.82 (26) 

0.52 

0.111±0.036 

4.1 

5 

142.733524 

E, 13i^i3-12i^i2 

49.3 

4.28(-5) 

2.06x1.67 (26) 

0.43 

0.190±0.042 

7.5 

2.06x1.67 (26) 

1.13 

0.195±0.044 

7.1 


142.735139 

A, 13i,i3-12i^i2 

49.3 

4.28(-5) 









6 

142.815476 

E, 13o,13-12o,12 

49.3 

4.28(-5) 

2.06x1.67 (26) 

0.30 

0.211±0.047 

7.6 

2.06x1.67 (26) 

1.20x0.70 (32) 

0.213±0.050 

8.3 


142.817021 

A, 13o,13-12o,12 

49.3 

4.28(-5) 









7 

142.924506 

E, 13i^i3-12o,12 

49.3 

6.59(-6) 

2.06x1.67 (26) 

point 

0.045±0.012 

8.8 

2.06x1.67 (26) 

point 

0.112±0.027 

9.7 


142.925911 

A, 13i,13-12o,12 

49.3 

6.59(-6) 









8 

164.955703 

E, 132,11-122,10 

58.5 

6.46(-5) 

2.31x1.72 (110) 

0.85 

0.159±0.051 

6.2 

2.39x1.77 (114) 

1.88 

0.236±0.058 

5.5 

9 

164.968638 

A, 132,11-122,10 

58.5 

6.46(-5) 

2.31x1.72 (110) 

0.4 

0.124±0.039 

3.5 

2.39x1.77 (114) 

0.38 

1 0.209±0.047 

7.0 

10 

163.829677 

E, 14i,13-13i,12 

62.5 

6.37(-5) 

2.31x1.72 (110) 

1.13x0.46 ( 66 ) 

0.131±0.056 

6.7 

2.28x1.76 

0.88 ^ 0.152±0.048 

4.7 

11 

163.835525 

A, 14i,i3-13i,12 

62.5 

6.37(-5) 

2.31x1.72 (110) 

1.82x1.36 (14) 

0.129±0.039 

3.6 

2.28x1.76 (-70) 

1.64x0.73 0.149±0.043 

4.4 

12 

165.653657 

E, 142,13-13i,12 

62.6 

7.63(-6) 

2.54x1.71 (113) 

1.6 

0.057±0.017 

6.3 

2.38x1.76 (114) 

1.65 

0.196±0.026 

14.3 


165.657529 

A, 142,13-13i,12 

62.6 

7.63(-6) 









13 

163.925845 

E, 15o,15-14i,14 

64.5 

1.02(-5) 

2.31x1.73 

1.89x1.30 (28) 

0.070±0.033 

7.8 

2.39x1.77 (114) (-61) 

2.40x0.86 

0.141±0.035 

6.2 


163.927369 

A, 15o,15-14i,14 

64.5 

1.02(-5) 









14 

163.960387 

A, 15i,15-14i^i4 

64.5 

6.53(-5) 

2.31x1.72 (110) 

1.48x0.56 (79) 

0.301±0.072 

6.9 

2.28x1.76 

0.96 

0.329±0.080 

5.9 


163.961884 

E, 15i,15-14i^i4 

64.5 

6.53(-5) 









15 

163.987455 

E, 15o,15-14o,14 

64.5 

6.54(-5) 

2.31x1.72 (110) 

0.78x0.33 (19) 

0.299±0.072 

6.9 

2.22x1.94 

0.93 

0.086±0.027 

6.0 

15 

163.988912 

A, 15o,15-14o,14 

64.5 

6.54(-5) 





2.22x1.94 

0.55 

0.136±0.034 


16 

164.022026 

E, 15i,i5-14o,14 

64.5 

1.03(-5) 

2.31x1.72 (110) 

1.81x0.76 (90) 

0.118±0.046 

9.7 

2.39x1.77 (114) 

point 

0.202±0.065 

10.0 


164.023416 

A, 15i,15-14o,14 

64.5 

1.03(-5) 









17 

164.205978 

E, 134 , 9 - 124^8 

64.9 

5.98(-5) 

2.31x1.72 (110) 

1.17 

0.126±0.039 

5.9 

2.28x1.76 

0.65 

0.193±0.088 

6.6 

18 

164.223815 

A, 134 , 9 - 124,8 

64.9 

5.98(-5) 

2.31x1.72 (110) 

1.29x0.90 (73) 

0.138±0.041 

6.5 

2.28x1.76 (0) 

3.02x1.74 

0.202±0.074 

5.2 

19 

142.664676 

E, 12i,ii-lli,io,'L’t —1 

234.0 

4.19(-5) 

2.06x1.67 (26) 

- 

< 0.040 

- 

2.06x1.67 (26) 

2.65x2.24 (0) 

0.037±0.009 

6.1 

20 

142.125411 

E, 13o,13-12o,125 —1 

236.9 

4.24(-5) 

2.06x1.67 (26) 

1.19 

0.033±0.011 

6.0 

2.06x1.67 (26) 

- 

< 0.048 

- 

21 

142.052774 

A, 13o,13-12o, 12 ) ^^t = l 

236.9 

4.21(-5) 

2.06x1.67 (26) 

- 

< 0.040 

- 

2.06x1.67 (26) 

2.24x0.56 (-29) 

0.038±0.018 

8.1 

22 

166.388878 

133 , 10 - 123 , 9 , vt = l 

248.0 

6.51(-5) 

2.54x1.71 (113) 

point 

0.072±0.032 

5.1 

2.63x1.73 

- 

< 0.043 

- 
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Table 11. Same as Table [Hlbut for CH 3 CN. 


IRAS2A IRAS4A 


N 

Frequency 

(GHz) 

Transition 

Eup 

(K) 

A„, 

(s-b 

Beam size 

("x", °) 

Source size 

("x". °) 

Flux 

(Jy km/s) 

dVw 

(km/s) 

Beam size 

("x", '’) 

Source size 

("x", °) 

Flux 

(Jy km/s) 

dVi\r 

(km/s) 

1 

165.569082 

9o-8o 

39.7 

2.65(-4) 

2.54x1.71 (113) 

0.99x0.78 (45) 

l.S45±0.373 

10.2 

2.38x1.76 (114) 

1.68x1.27 (-67) 

1.058±0.217 

1 9.4 


165.565891 

9i-8i 

46.9 

2.62(-4) 








Cn 

2 

165.556322 

92-82 

68.3 

2.52(-4) 

2.54x1.71 (113) 

1.07x0.70 (-57) 

0.S1S±0.170 

5.5 

2.38x1.76 (114) 

0.55 

0.533T0.116 

° 5.1 

3 

165.540377 

93-83 

104.0 

2.36(-4) 

2.54x1.71 (113) 

0.98x0.77 (65) 

0.986±0.205 

6.1 

2.38x1.76 (114) 

0.51 

0.493±0.104 

' 5.7 

4 

165.518064 

94-84 

154.0 

2.13(-4) 

2.54x1.71 (113) 

0.87x0.79 (-26) 

0.541±0.116 

6.1 

2.38x1.76 (114) 

1.22x0.5 (-56) 

0.371±0.081 

6.0 

5 

165.489391 

95-85 

218.3 

1.83{-4) 

2.54x1.71 (113) 

0.44 

0.426±0.097 

6.8 

2.38x1.76 (114) 

0.52 

0.482T0.109 

8.7 

6 

165.454370 

95-85 

296.8 

1.47(-4) 

2.54x1.71 (113) 

0.74x0.45 (-80) 

0.331±0.078 

6.6 

2.38x1.76 (114) 

1.29x0.65 (-74) 

0.310T0.070 

5.3 

7 

165.413015 

97-87 

389.5 

1.04(-4) 

2.54x1.70 

0.44 

0.111±0.049 

6.4 

2.38x1.76 (114) 

1.88x0.75 (-20) 

0.075±0.032 

4.8 








Table 12. Same as Table [8] but for CH3OCH3 


N 

Frequency 

(GHz) 

Transition 

Eup 

(K) 

A„i 

(s“b 


IRAS2A 




IRAS4A 



Beam size 

("x", °) 

Source size 

("x".°) 

Flux 

(Jy km/s) 

(km/s) 

Beam size 

("x", “) 

Source size 

("x", °) 

Flux 

(Jy km/s) 

dViY 

(km/s) 

1 

143.017994 

32,2"2i,i , EA 

11.1 

1.09e-5 

2.06x1.67 (26) 

0.81 

0.044±0.023 

10.1 

2.06x1.67 (26) 

3.81x1.66 (89) 

0.060±0.022 

8.5 


143.018373 

32,2-2i,i, AE 

11.1 

1.09e-5 










143.020764 

32,2"2i,i , EE 

11.1 

1.09e-5 










143.023345 

32,2-2i,i, AA 

11.1 

1.09e-5 









2 

144.856766 

63,3-62,4 ea 

31.8 

1.04(-5) 

2.08x1.65 (30) 

0.38 

0.070±0.016 

5.1 

2.16x1.73 (25) 

1.29x1.28 (0) 

0.161T0.045 

19.6 


144.858984 

63,3-62,4 ee 

31.8 

1.06(-5) 










144.855091 

63,3-62,4 ae 

31.8 

1.06(-5) 









3 

143.599420 

"^3,4-72,5 AE 

38.2 

1.10(-5) 

2.08x1.65 (30) 

1.39x0.94 (-11) 

0.083±0.030 

16.1 

2.16x1.73 (25) 

point 

0.136±0.044 

14.0 


143.600084 

63,3-62,4 ea 

38.2 

1.10(-5) 










143.602993 

63,3-62,4 ee 

38.2 

1.10(-5) 










143.606232 

63,3-62,4 aa 

38.2 

1.10(-5) 









4 

141.828855 

33,5-82,6! AE 

45.5 

l.lle-5 

2.06x1.67 (26) 

0.65 

0.149±0.045 

14.3 

2.06x1.67 (26) 

0.46 

0.140±0.049 

16.7 


141.829146 

83,5-82,6 ! ea 

45.5 

l.lle-5 







1 



141.832261 

83,5-82,6 ! ee 

45.5 

l.lle-5 







51 



141.835521 

83,5-82,6 ! aa 

45.5 

l.lle-5 







1 


5 

143.159951 

132,12-13i^13, ea 

88 

7.75e-6 

2.06x1.67 (26) 

0.31 

0.050±0.012 

7.1 

2.06x1.67 (26) 

point 

0.079±0.022 

13.0 


143.159952 

132,12-13i^13, AE 

88 

7.75e-6 










143.162986 

132,12-13i^13, ee 

88 

7.75e-6 










143.166020 

132,12-13i,13, aa 

88 

7.75e-6 









6 

165.208844 

163,13-152,14 ea 

122 

1.75(-5) 

2.31x1.72 (110) 

point 

0.121±0.046 

13.8 

2.39x1.77 (114) 

2.47x1.20 (-26) 

0.115±0.025 

9.5 


165.208848 

163,13-152,14 ae 

122 

1.75(-5) 










165.211731 

163,13-162,14 ee 

122 

1.75(-5) 










165.214617 

163,13-162,14 aa 

122 

1.75(-5) 









7 

164.988708 

203,18-194,15 AA 

204 

4.60(-6) 

2.31x1.72 (110) 

point 

0.074±0.033 

5.5 

2.39x1.77 (114) 

- 

< 0.072 

- 


164.990831 

203,18-194,15 EE 

204 

4.60(-6) 










164.992951 

203,18-194,15 EA 

204 

4.60(-6) 










164.992958 

203,18-194,15 AE 

204 

4.60(-6) 









8 

142.403201 

25i,24-252,23 ) ea 

313.5 

1.57e-5 

2.06x1.67 (26) 

point 

0.051±0.013 

7.2 

2.06x1.67 (26) 

0.41 

0.035±0.011 

6.0 


142.403201 

25i,24-252,23' 

AE 

313.5 

1.57e-5 

142.404442 

251,24-252,23' 

EE 

313.5 

1.57e-5 

142.405682 

251,24-252,23' 

AA 

313.5 

1.57e-5 
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N 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

13 

14 


Table 13. Same as Table [ 8 ] but for C 2 H 5 OH. 


IRAS2A IRAS4A 


Frequency 

Transition 

Eiip 

^ul 

Beam size 

Source size 

Flux 

dVw 

Beam size 

Source size 

Flux 

dV\\r 

(GHz) 


(K) 

(s“b 

("x", °) 


(Jy km/s) 

(km/s) 

("x", °) 

("x". °) 

(Jy km/s) 

(km/s) 

142.285054 

9o, 9 , 2 - 81 , 8,2 

37.2 

1.51{-5) 

2.06x1.67 (26) 

2.74x0.47 (1) 

0.027±0.011 

4.1 

2.06x1.67 (26) 

0.52 

0.045±0.015 

8.2 

164.900973 

60 , 6 , 4,0 

78.8 

1.22(-5) 

2.31x1.72 (110) 

- 

< 0.066 

- 

2.39x1.77 (114) 

- 

< 0.072 

- 

144.057496 

133,11,2-132,12,2 

87.9 

1.80(-5) 

2.08x1.65 (30) 

point 

0.031±0.010 

- 

2.16x1.73 (25) 

point 

0.067T0.018 

4.1 

141.820317 

81 ,7,0-’^!, 6,0 

88.8 

2.46e-5 

2.06x1.67 (26) 

1.16 

0.034±0.013 

4.1 

2.06x1.67 (26) 

point 

0.039±0.015 

6.4 

164.626167 

54 , 1 , 0 - 43 , 1,1 

88.8 

2.06(-5) 

2.31x1.72 (110) 

- 

< 0.066 

- 

2.39x1.77 (114) 

- 

< 0.072 

- 

164.630894 

54 , 2 , 0 - 43 , 2,1 

88.8 

2.06(-5) 

2.31x1.72 (110) 

- 

< 0.066 

- 

2.39x1.77 (114) 

- 

< 0.081 

53 

144.493107 

142,13,2-1411,4,2 

92.6 

1.86(-5) 

2.08x1.65 (30) 

- 

< 0.046 

- 

2.16x1.73 (25) 

1.19 

0.045±0.021 

1 

166.259891 

10i,io,o-9i,9,o 

102.1 

4.03(-5) 

2.54x1.71 (113) 

point 

0.114±0.049 

7.3 

2.38x1.76 (114) 

1.08 

O.II 81 I 1 O.O 53 

4.2 

166.758214 

10i,io,i-9i,9,i 

106.8 

4.07(-5) 

2.54x1.71 (113) 

0.95 

0.043±0.030 

- 

2.38x1.76 (114) 

point 

0.053±0.020 

- 

142.083012 

ll2,10,l-lll,10,0 

121 

9.96(-6) 

2.06x1.67 (26) 

point 

0.028±0.009 

16.3 

2.06x1.67 (26) 

- 

< 0.040 

- 

164.511879 

183 , 15 , 2 - 174 , 14,2 

156.8 

8.73(-6) 

2.31x1.72 (110) 

point 

0.036±0.021 

3.6 

2.39x1.77 (114) 

- 

< 0.072 

- 

164.429108 

143 , 12 , 1 - 134 , 10,0 

160.0 

4.25(-6) 

2.31x1.72 (110) 

- 

< 0.066 

- 

2.39x1.77 (114) 

- 

< 0.072 

- 

142.046310 

202,18,2-20i,19,2 

185.5 

1.77(-5) 

2.06x1.67 (26) 

- 

< 0.040 

- 

2.06x1.67 (26) 

- 

< 0.040 

- 

163.515766 

213 , 19 , 2 - 204 , 16,2 

205.4 

9.22(-6) 

2.31x1.72 (110) 

- 

< 0.066 

- 

2.39x1.77 (114) 

- 

< 0.072 

- 

141.735490 

214 , 17 , 2 - 213 , 18,2 

215.5 

1.72e-5 

2.06x1.67 (26) 

point 

0.030±0.015 

8.0 

2.06x1.67 (26) 

point 

0.048±0.018 

7.6 








Table 14. Same as Table [H] but for HCOCH 2 OH. 


IRAS2A IRAS4A 


N 

Frequency 

Transition 

Eup 

A„i 

Beam size 

Source size 

Flux 

dVy[r 

Beam size 

Source size 

Flux 

dVy/ 


(GHz) 


(K) 

(b-M 

("x", °) 

("x". °) 

(Jy km-s) 

(km/s) 

("x",°) 

("x", °) 

(Jy km-s) 

(km/s) 

1 

163.951686 

83 , 5 - 72 , 6 ! 1^ — 0 

25.6 

3.68(-5) 

2.31x1.72 (110) 

- 

< 0.066 

- 

2.39x1.77 (114) 


- 

- 

2 

143.640947 

14o,14-13i,13, = 0 

53.1 

8.04(-5) 

2.08x1.65 (30) 

point 

0.037±0.020 

- 

2.06x1.67 (26) 

- 

< 0.129 

- 

3 

143.765755 

14i,14-13o,13i = 0 

53.1 

8.06(-5) 

2.08x1.65 (30) 

- 

< 0.046 

- 

2.06x1.67 (26) 

point 

0.069±0.035 

S 5.6 

4 

164.047038 

152,14-14i,i 3, r) = 0 

66.5 

9.25(-5) 

2.31x1.72 (110) 

point 

0.076±0.034 

7.4 

2.39x1.77 (114) 

1.89 

0.132±0.051 

1 4.5 

5 

163.542260 

16o,16-15i, 15 , r) = 0 

68.3 

1.21(-4) 

2.31x1.72 (110) 

point 

0.037±0.023 

7.8 

2.39x1.77 (114) 

- 

0.172±0.050 

8.0 

6 

163.580057 

16i,16-15o,15i = 0 

68.3 

1.21(-4) 

2.31x1.72 (110) 

point 

0.038±0.021 

8.5 

2.39x1.77 (114) 

point 

0.103±0.038 

5.0 

7 

163.697251 

127 , 5 - 126 ^ 6 , 'u — 0 

73.1 

5.60(-5) 

2.31x1.72 (110) 

- 

< 0.066 

- 

2.39x1.77 (114) 

point 

0.071±0.028 

3.6 

8 

163.709163 

127 ,6-126,7i V = 0 

73.1 

5.60(-5) 

2.31x1.72 (110) 

- 

< 0.066 

- 

2.39x1.77 (114) 

1.55 

0.076±0.045 

3.8 

9 

142.784665 

236,18-235,19 

177 

5.19(-5) 

2.08x1.65 (30) 

- 

< 0.040 

- 

2.06x1.67 (26) 

point 

0.039±0.014 

5.2 








N 

1 

2 

3 

4 

5 

6 

8 

7 

9 

10 


Table 15. Same as Table [ 8 ] but for C 2 H 5 CN. 


IRAS2A IRAS4A 


Frequency 

Transition 

E^p 

A„, 

Beam size 

Source size 

Flux 

dV\y 

Beam size 

Source size 

Flux 

dV^Wr 

(GHz) 


(K) 

(s"b 

("x", °) 

("x", “) 

(Jy km/s) 

(km/s) 

("x", °) 

("x", ») 

(Jy km/s) 

(km/s) 

142.346330 

162,15-152,14 

62.7 

2.37(-4) 

2.06x1.67 (26) 

2.33 

0.029±0.010 

12.0 

2.06x1.67 (26) 

2.51x1.17 (-40) 

0.045±0.013 

6.8 

143.529200 

163 , 14 - 153,13 

68.5 

2.39(-4) 

2.08x1.65 (30) 

point 

0.025±0.009 

4.1 

2.16x1.73 (25) 

- 

< 0.081 

- 

144.104740 

163 , 13 - 153,12 

68.6 

2.42(-4) 

2.08x1.65 (30) 

- 

< 0.046 

- 

2.16x1.73 (25) 

- 

< 0.081 

- 

143.506970 

164 , 13 - 154,12 

76.3 

2.32(-4) 

2.08x1.65 (30) 

- 

< 0.046 

- 

2.16x1.73 (25) 

- 

< 0.081 

- 

143.535290 

164 , 12 - 154,11 

76.3 

2.32(-4) 

2.08x1.65 (30) 

- 

< 0.046 

- 

2.16x1.73 (25) 

- 

< 0.081 

1 - 

164.584755 

19o,19-18o,18 

80 

3.74 (-4) 

2.31x1.72 (110) 

- 

< 0.066 

- 

2.39x1.77 (114) 

- 

< 0.072 

55 

163.948705 

CD 

00 

80.1 

3.69(-4) 

2.31x1.72 (110) 

- 

< 0.066 

- 

2.39x1.77 (114) 

point 

0.110ih0.044 

1 8.5 

143.406554 

165 , 12 - 155,11 

86.3 

2.23(-4) 

2.08x1.65 (30) 

point 

0.024±0.016 

5.8 

2.16x1.73 (25) 

- 

< 0.081 

- 

143.407188 

165,11-155,10 

86.3 

2.23(-4) 

2.08x1.65 (30) 







- 

143.335284 

168,8-158,7 

129.6 

1.85(-4) 

2.06x1.67 (26) 

point 

0.045±0.011 

7.4 

2.06x1.67 (26) 

0.17 

0.037±0.012 

11.2 

143.335284 

168,9-158,8 

129.6 

1.85(-4) 









143.337710 

167,10-157,9 

112.9 

1.99(-4) 









143.337710 

167,9-157,8 

112.9 

1.99(-4) 









143.343925 

169,7-159,6 

148.4 

1.68(-4) 

2.06x1.67 (26) 

point 

0.034±0.010 

7.1 

2.06x1.67 (26) 

1.43 

0.045±0.012 

4.1 

143.343925 

169,8-159,7 

148.4 

1.68(-4) 
















Table 16. Same as Table [H] but for other molecules. 


IRAS2A IRAS4A 


Molecule 

Frequency 

Transition 

Eup 

A„i 

Beam size 

Source size 

Flux 

dVyy 

Beam size 

Source size 

Flux 

dVy\/ 


(GHz) 


(K) 

(b- 1 ) 

("x", °) 

("x". °) 

(Jy km/s) 

(km/s) 

("x", °) 

("x", °) 

(Jy lim/s) 

(km/s) 

HC 3 N 

163.753389 

18-17 

74.7 

3.45(-4) 

2.31x1.72 (110) 

2.66x1.80 (34) 

0.572±0.116 

6.2 

2.39x1.77 (114) 

6.03x1.96 (17) 

0.55^D.121 

7.4 


141.98374 

2 o, 2 - 1 o ,1 

10.2 

7.25(-5) 

2.06x1.67 (26) 

0.60 

0.115±0.02S 

6.3 

2.06x1.67 (26) 

2.45x1.25 (0) 

0.062^0.022 

7.3 


143.213062 

2i,i-li,o 

22.2 

5.57(-5) 

2.06x1.67 (26) 

1.94x1.05 (1) 

0.028±0.007 

- 

- 

- 

- 

- 

NH 2 CHO 

142.701479 

77 , 1 , 7 , 8 - 61 , 6,7 

30.4 

2.02(-4) 

2.06x1.67 (26) 

point 

0.185±0.043 

8.2 

2.06x1.67 (26) 

point 

0.090±0.025 

7.4 

CH 2 CO 

142.76892 

7i,6-6i,5 

40.5 

3.10(-5) 

2.06x1.67 (26) 

1.03x0.61 (30) 

0.119±0.028 

5.8 

2.06x1.67 (26) 

1.80x1.34 (34) 

0.184±0.040 

5.1 
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Fig. 9.— Rotational and population diagrams of methanol isotopologues (CH 3 OH, 
^^CHsOH), HCOOCH 3 , and CH 3 CN for source sizes derived from the PD analysis of the 
methanol population distribution (0.36 " for IRAS2A and 0.20 " for IRAS4A). Observa¬ 
tional data is depicted by the black diamonds. Error bars are derived assuming a calibration 
uncertainty of 20 % on top of the statistical error. Straight lines represent the best £t of the 
RD analysis to the data. Red crosses show the best fit of the PD to the data. 
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Fig. 10.— Rotational and population diagrams of CH 3 OCH 3 , C 2 H 5 OH, HCOCH 2 OH, 
C 2 H 5 CN for source sizes derived from the PD analysis of the methanol population dis¬ 
tribution (0.36 " for IRAS2A and 0.20 " for IRAS4A). Observational data is depicted by the 
black diamonds. Error bars are derived assuming a calibration uncertainty of 20 % on top 
of the statistical error. Straight lines represent the best £t of the RD analysis to the data. 
Red crosses show the best £t of the PD to the data. 






































Table 17. COMs abundances observed in high-mass, intermediate-mass, and low-mass hot 

cores. 





CH 3 OH 

H 2 CO 

CH 2 CO 

HCOOCH 3 

HCOCH 2 OH 

CH 3 OCH 3 

C 2 H 5 OH 

CH 3 CN 

C 2 H 5 CN 

Source 

^bol 

JV(H 2 ) 

Trot 

^H2 

-^meth 

■^meth 

-^meth 

-^meth 

-^meth 

■^meth 

2 ^meth 

2 ^meth 


(Lq) 

(cm- 2 ) 

(K) 














High-mass protostars observed with 

single-dish telescopes 




IRAS 20126+4104 

1.3(+4) 

1.0( + 24) 

300 

l.l(-7) 

2.9(-l) 

< l.l(-2) 

< 1.8(-2) 

- 

< 9.1(-2) 

< 2.7(-2) 

5.2(-2) 

< 4.5(-3) 

IRAS 18089-1732 

3.2(+4) 

1.0( + 24) 

300 

2.0(-7) 

8.6{-2) 

3.1(-2) 

1.3(-1) 

- 

5.9(-l) 

1.8(-1) 

2.1(-2) 

1.4(-2) 

G31.41+0.31 

2.6(+5) 

1.7(+23) 

200 

6.0(-6) 

5.4{-2) 

2.8(-2) 

1.4(-1) 

i.i(-i) 

4.5(-l) 

1.3(-1) 

1.3(-2) 

1.7(-2) 

AFGL 2591 

2.0(+4) 

7.6(+22) 

147 

6.2(-7) 

2.8(-l) 

2.1(-2) 

< 5.1(-1) 

- 

< 1.6(-1) 

< 2.1(-2) 

< 7.4(-2) 

< 1.6(-2) 

NGC 7538 IRSl 

1.3(+5) 

2.1(+23) 

156 

5.7(-7) 

2.1(-1) 

5.3(-2) 

1.2(-1) 

- 

< 1.3(-1) 

4.8(-2) 

< 6.8(-2) 

< 7.7(-3) 

G24.78 

7.9(+5) 

4.0(+23) 

211 

7.0(-7) 

2.3{-l) 

5.2(-2) 

l.l(-l) 

- 

4.3(-l) 

2.5(-2) 

2.1(-1) 

1.4(-2) 

G75.78 

1.9(+5) 

1.2(+23) 

113 

9.2(-7) 

2.0{-l) 

6.2(-2) 

6.5(-2) 

- 

2.1(-1) 

< 2.2(-2) 

1.6(-2) 

< 

W33A 

1.0(+5) 

2.6(+23) 

259 

7.7(-7) 

2.7{-l) 

4.9(-2) 

1.3(-1) 

- 

1.4(-1) 

2.4(-2) 

1.4(-1) 

< l-l(-2) 

NGC 6334 IRSl 

1.7(+5) 

2.4(+23) 

178 

4.0(-6) 

1.3(-1) 

2.0(-2) 

1.2(-1) 

- 

6.0(-l) 

2.0(-2) 

3.0(-2) 

5.3(-3) 

W3 (H20) 

2.0(+4) 

1.8(+23) 

181 

5.6(-6) 

1.8(-1) 

1.5(-2) 

5.2(-2) 

- 

1.5(-1) 

8.4(-3) 

7.0(-3) 

4.5(-3) 

Sgr B2 (M) 

6.5(+6) 

3.5(+24) 

150 

7.4(-9) 

l.O(-l) 

1.9(-2) 

6.2(-2) 

- 

- 

4.2(-2) 

6.0(-2) 

- 

Sgr B2 (N) 

6.5(+6) 

8.0(+24) 

170 

6.3(-7) 

4.8(-2) 

1.4(-2) 

< 3.8(-2) 

- 

5.6(-2) 

3.8(-2) 

5.1(-2) 

1.3(-1) 

G327.3-0.6 

1.0(+5) 

3.0(+24) 

118 

2.0(-5) 

3.6(-5) 

- 

8.0(-2) 

- 

5.4(-l) 

4.1(-3) 

3.5(-2) 

2.2(-2) 

Orion KL - HC 

1.0(+5) 

3.1(+23) 

128 

2.2(-6) 

5.5(-2) 

- 

- 

- 

3.1(-2) 

- 

1.4(-2) 

5.0(-3) 

Orion KL - CR 

1.0(+5) 

3.9(+23) 

140 

1.2(-6) 

3.7{-2) 

4.3(-3) 

2.8(-l) 

- 

1.4(-1) 

1.4(-2) 

l.l(-2) 

- 

G34.3+0.15 

6.3{+5) 

5.3(+23) 

336 

7.0(-8) 

- 

1.8{-2) 

4.3(-l) 

- 

2.5(-l) 

9.5(-2) 

6.5(-3) 

- 

G34.3+0.2 

6.3(+5) 

1.6(+23) 

96 

1.7(-7) 

- 

- 

5.4(-2) 

- 

1.4(-1) 

6.5(-2) 

- 

1.0(-2) 

DR21(OH) 

5.0(+4) 

2.5(+24) 

150 

1.0(-8) 

- 

- 

< 4.0(-3) 

- 

- 

< 5.2(-3) 

- 

< 8.4(-4) 

W51 

1.5(+6) 

3.3(+23) 

208 

3.0(-7) 

- 

- 

1.2(-1) 

- 

- 

3.1(-2) 

- 

7.0(-3) 






High-mass 

hot cores 

observed with interferometers 




Orion KL 

1.0(+5) 

4.4(+24) 

200 

4.5(-9) 

- 

- 

3.0(-l) 

- 

5.0(-l) 

l.O(-l) 

l.O(-l) 

2.5(-l) 

G29.96 

9.0(+4) 

3.3(+24) 

200 

1.2(-7) 

- 

- 

2.0(-l) 

- 

5.0(-l) 

1.5(-1) 

2.5(-2) 

2.5(-2) 

G19.61-0.23 

1.6(+5) 

8.4(+23) 

151 

6.2(-7) 

- 

- 

4.2(-2) 

- 

2.7(-2) 

1.2(-1) 

7.9(-2) 

3.1(-2) 





Intermediate-mass hot cores observed with interferometers 




I22198-MM2 

370 

2.0(+25) 

120 

1.2(-6) 

- 

- 

l.l(-2) 

- 

- 

9.1(-3) 

- 

- 

A5142-MM1 

2300 

1.0(+25) 

210 

2.3(-7) 

- 

- 

< 8.7(-2) 

- 

- 

9.1(-2) 

- 

- 

A5142-MM2 

2300 

2.0(+25) 

140 

2.0(-7) 

- 

- 

< 5.0(-2) 

- 

- 

2.0(-2) 

- 

- 

NGC 7129 FIRS2 

500 

2.5(+24) 

238 

1.0(-6) 

1.6(-2) 

- 

1.5(-2) 

- 

1.2(-2) 

8.8(-3) 

5.2(-3) 

3.5(-4) 





Low-mass protostars observed with 

single-dish telescopes 




IRAS 16293 

27 

2.0(+23) 

84 

1.0(-7) 

2.6(-l) 

1.0(-3) 

9.0(-2) 

6.9(-3) 

4.0(-l) 

< 5.0(-2) 

9.1(-3) 

< 2.0(-3) 

IRAS2A 

36 

2.1(+23) 

101 

8.8(-7) 

2.9{-l) 

- 

< 8.5(-l) 

- 

< 5.3(-l) 

- 

l.l(-2) 

< 1.3(-1) 

IRAS4A 

9.1 

1.6(+24) 

24 

1.4(-7) 

1.4(-1) 

- 

5.5(-l) 

- 

< 2.2(-l) 

- 

1.3(-2) 

< 9.2(-3) 

IRAS4B 

4.4 

8.1(+22) 

34 

6.9(-6) 

2.0(-l) 

- 

1.3(-1) 

- 

< 1.9(-1) 

- 

1.6(-2) 

< 1.2(-1) 

SMMl 

30 

1.3(+23) 

16 

1.9(-9) 

- 

- 

l.O(-l) 

- 

5.3(-2) 

< 3.4(-2) 

- 

- 


Reference 


1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

3 

4 



9 

10 
11 

12 

12 

12 

13 

14 

15 
15 

15 

16 



Table 17—Continued 


SMM4 

1.9 

l.l(+23) 

13 

9.5(-9) 

- 

- 

< 1 . 0 (- 2 ) 

- 

< 8.0(-3) 

< 6.0(-3) 

- 

- 

16 

Bl-a 

1.3 

1.9(-b22) 

15 

6.4(-10) 

- 

- 

l.O(-l) 

- 

< 6.7(-2) 

- 

1.3(-2) 

- 

17 

SVS 4-5 

38 

5.7(-b22) 

20 

3.9(-9) 

- 

- 

4.5(-2) 

- 

l.O(-l) 

- 

7.7(-3) 

- 

17 

B5 IRSl 

4.7 

2.3(-b22) 

17 

1.0(-9) 

- 

- 

< 1.7(-1) 

- 

< 3.5(-l) 

- 

1.7(-2) 

- 

17 

IRAS 03235 

1.9 

1.4(+23) 

18 

8 . 6 (-ll) 

- 

- 

< 1.7(-1) 

- 

< 4.2(-l) 

- 

< 2.5(-2) 

- 

17 

IRAS 04108 

0.62 

2.9(-b22) 

9 

4.1(-10) 

- 

- 

< 8.3(-2) 

- 

- 

- 

- 

- 

17 

L1489 IRS 

3.7 

4.3(-b22) 

8 

1 . 2 (- 10 ) 

- 

- 

< 4.0(-l) 

- 

- 

- 

- 

- 

17 





Low-mass hot 

corinos observed with 

interferometers 





IRAS2A-RD 

36 

2.0(-b24) 

179 

2.6(-7) 

8 . 1 (- 2 ) 

1.4(-3) 

1.9(-2) 

1.5(-3) 

1 . 2 (- 2 ) 

1.5{-2) 

3.0(-3) 

2.7(-4) 

18 

IRAS2A-PD 

36 

2.0(-b24) 

140 

1 . 0 (- 6 ) 

- 

- 

1 . 6 (- 2 ) 

1.4(-3) 

1 . 0 (- 2 ) 

1 . 6 {- 2 ) 

4.0(-3) 

3.0(-4) 

18 

IRAS4A-RD 

9.1 

1.4(+25) 

300 

1.7(-8) 

6.3(-2) 

2.8(-3) 

1.5(-2) 

2.5(-3) 

8.7(-3) 

1 . 2 (- 2 ) 

1.8(-3) 

4.2(-4) 

18 

IRAS4A-PD 

9.1 

1.4(+25) 

140 

4.3(-7) 

- 

- 

3.1(-2) 

3.0(-3) 

1 . 0 (- 2 ) 

1 . 0 (- 2 ) 

3.9(-3) 

4.0(-4) 

18 


Note. — 1: 

Isokoski et al. 1 2013l'l: 2: IpisschoD et al.l 120071): 3: 

Nummelin et al.l 1 

200d'): 4: iNeill et al. i2014l'): 5: Icibb et al.N200cl^: 6: 

Crockett et al.li2014h 7: MacDonald et ai. il996h: 8 

Ikeda et al. 

2001 

): 9: Beuther et al 

l200fl 

): 10: Beuther et a 

12007 

): 11: Oin et al. 

boiof); 

12: Palau et al. 

120111): 13:lFuente et al. 120141^: 14: 

Maret et al. 

20051^ :|jaber et al.l i201^ 

); 15: 

Maret et al.l i2004l. 

200,Jl: 

Bottinelli et al.l 

J2004al. 

2 OO 7 I 1 : 16: Obere- et al.l J 2 OI 1 I 1 : 17: Oberp; et al.| J2014|'t: 18: This work 



















































